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Fractures  in  a  rock  mass  can  be  considered  as  pre-existing  flaws  in  an  otherwise 
continuous  medium.  It  was  initially  thought  that  crack  propagation  criteria  from 
fracture  mechanics,  which  were  experimentally  verified  on  polymers,  could  be 
applied  to  this  conceptual  model  of  fractured  rock.  However,  a  systematic 
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the  basis  of  fracture  mechanics  theory  and  what  had  been  shown  in  experiments 
on  polymers.  The  most  significant  observation  was  the  coalescence  of  pre¬ 
existing  flaws  through  secondary  crack  growth  that  occurred  after  the  usual 
appearance  of  wing  cracks. 

Linear  elastic  analyses  of  the  experimental  flaw  geometries  indicated  that  a 
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Chapter  1 
Introduction 


Rock  mass  failure  has  been  viewed  as  occuring  in  either  one  of  two  ways:  (1) 
failure  along  continuous  pre-existing  fractures,  or  (2)  fracturing  of  the  intact  rock.  In 
frequent  cases  where  the  fractures  are  not  continuous  (Fig.  1.1a),  failure  occurs  through 
fracturing  of  the  intact  rock  bridges  that  leads  to  coalescence  of  the  pre-existing  fractures 
(Fig.  1.1b).  Present  models  for  rock  mass  failure  prediction  typically  assume  that  the 
fractures  are  continuous,  thus  neglecting  the  significant  contribution  of  the  intact  rock 
bridges  to  the  overall  strength  of  the  rock  mass.  The  objective  of  this  study  is  to  develop  a 
mechanical  model  for  fracture  coalescence  that  can  be  used  to  improve  failure  prediction  of 
rock  masses  with  discontinuous  fractures. 

Fractures  in  a  rock  mass  can  be  considered  as  pre-existing  cracks  in  an  otherwise 
continuous  medium.  It  was  initially  thought  that  crack  propagation  criteria  from  fracture 
mechanics,  which  have  been  verified  using  polymeric  materials,  could  be  applied  to  this 
conceptual  model  of  fractured  rock  (Chan,  1986).  However,  published  results  of 
experiments  exhibited  crack  growth  in  rock  that  was  very  different  from  what  was  seen  in 
polymers  (see  Chapter  2).  Preliminary  tests  using  gypsum  as  a  model  material  for  rock 
showed  coalescence  between  pre-existing  fractures  that  had  never  occured  in  similar  tests 
on  polymers.  More  importantly,  the  observed  coalescence  could  not  be  predicted  by 
existing  crack  propagation  criteria.  Since  gypsum  more  closely  modelled  rock  behavior, 
the  fracture  patterns  in  the  preliminary  tests  were  more  likely  to  occur  in  a  fractured  rock 
mass.  Thus,  these  initial  findings  led  to  a  systematic  experimental  and  analytical 
investigation  of  crack  growth  and  coalescence  mechanisms  in  rock,  using  gypsum  as  a 
model  material.  The  study  consisted  of: 
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(1)  Monitoring  crack  growth  in  uniaxially  compressed  gypsum  specimens  with  pre-existing 
fractures  with  various  arrangements.  Consistent  coalescence  behavior  as  a  function  of  pre¬ 
existing  fracture  geometry  was  noted  and  is  documented  in  Chapter  3. 

(2)  Development  of  an  analytical  coalescence  model  based  on  correlations  between 
computed  stress/strain  fields  and  observed  fracturing  in  the  experiments.  Using  this 
analytical  model,  satisfactory  agreement  was  achieved  between  simulated  and  experimental 
coalescence  behavior  as  a  function  of  pre-existing  fracture  geometry.  The  development, 
application  and  validation  of  this  coalescence  model  are  described  in  Chapter  4. 

Conclusions  and  recommendations  for  future  study  are  then  given  in  Chapter  5. 

Note:  In  the  following  chapters  of  this  report,  the  term  "flaw"  will  be  used  to  refer  to  pre¬ 
existing  fractures.  Also,  the  standard  mechanics  sign  convention  is  used,  i.e.,  positive  and 
negative  values  refer  to  tensile  and  compressive  stress/strain  components,  respectively. 
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Pre-existing  fractures 


Figure  1.1(a)  Rock  slope  with  discontinuous  fractures 


Figure  1.1(b)  Failure  in  a  rock  slope  due  to  coalescence  of  pre-existing  fractures 


Chapter  2 

Previous  Work:  Crack  Growth  from  Inclined  Flaws 
in  Compression 


2.1  Introductory  Remarks 

Starting  with  the  pioneering  work  of  Erdogan  and  Sih  (1963),  numerous  papers 
have  been  written  on  the  problem  of  predicting  crack  growth  from  inclined  pre-existing 
flaws,  also  referred  to  as  mixed-modf  crack  propagation.  Comprehensive  reviews  of  the 
various  crack  propagation  criteria  and  supporting  experimental  studies  have  already  been 
presented  by  Maiti  and  Smith  (1984)  and  by  Yarema  and  Ivanitskaya  (1986).  There  has 
also  been  an  entire  conference  devoted  to  mixed-mode  crack  propagation  (Sih  and 
Theocaris,  1981)  as  well  as  a  book  written  by  Gdoutous  (1984)  on  the  same  subject.  Since 
such  comprehensive  reviews  already  exist,  this  chapter  will  focus  on  work  reported  in  the 
literature  that  is  most  relevant  to  this  study,  i.e.,  crack  propagation  from  inclined  flaws 
subjected  to  far-field  compression  and,  most  importantly,  coalescence  between  multiple 
pre-existing  flaws. 

The  first  section  of  this  chapter  presents  experimental  work  performed  by  others  on 
(1)  crack  growth  from  single  flaws,  and  (2)  crack  coalescence  between  multiple  flaws. 
Two  general  types  of  materials  were  used  in  these  experiments:  polymers  such  as  glass, 
PMMA  and  CR39,  and  rocks.  Polymers  are  favored  model  materials  for  studying  rock 
structures  due  to  their  photoelastic  properties  and  easy  machinability.  However,  this 
review  will  show  that  conclusions  drawn  from  fracture  experiments  in  polymers  may  not 
be  applicable  to  fracturing  of  rock.  Specifically,  differences  between  crack  propagation 
from  single  flaws  in  polymers  and  in  rocks  can  be  clearly  established  by  comparing 
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published  results  between  several  experiments  using  these  materials.  On  the  other  hand,  as 
far  as  crack  coalescence  is  concerned,  only  one  such  experiment  is  known  to  have  been 
performed  in  rocks  and  its  results  differed  significantly  from  similar  tests  on  polymers. 
The  observations  in  this  single  experiment  corroborate  the  findings  presented  in  Chapter  3 
of  this  report. 

The  second  section  of  this  chapter  describes  analyses  using  existing  crack 
propagation  criteria  that  successfully  modelled  the  observed  crack  growth  from  flaws  in 
polymer  specimens.  However,  the  observed  crack  growth  from  flaws  in  rock  was  not 
explicitly  simulated  by  these  analyses,  thus  highlighting  the  inadequacy  of  existing  crack 
propagation  criteria  for  modelling  rock  fracture. 
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2.2  Observed  Crack  Growth  from  Single  and  Multiple  Flaws 


2.2.1  Specimens  with  Single  Inclined  Flaws 

It  has  been  been  well  documented  in  the  literature  that  when  a  single,  inclined  flaw 
is  subjected  to  uniaxial  compression,  cracks  initiate  from  the  flaw  tips  approximately 
perpendicular  to  the  plane  of  the  flaw,  curve  and  propagate  parallel  to  the  applied  uniaxial 
compressive  load  (see  Fig.  2.1).  These  cracks,  referred  to  as  wing  cracks  in  this  study, 
grow  stably,  i.e.,  an  increase  in  load  is  required  for  them  to  lengthen.  Wing  cracks  were 
seen  growing  from  inclined  flaws  in  polymer  materials  such  as  polymethymethacrylate  or 
PMMA  (McClintock,1963;  Adams  and  Sines,  1978;  Ingraffea  and  Heuze,  1980),  glass 
(Brace  and  Bombolakis,  1963;  Hoek  and  Bieniawski,  1965)  and  Columbia  Resin  39  or 
C7?J9(Horii  and  Nemat-Nasser,  1986).  Specimens  with  open  and  closed  pre-existing 
flaws  both  exhibited  wing  crack  growth. 

Wing  cracks  were  also  seen  in  uniaxial  compression  tests  on  rock  specimens  with 
pre-existing,  open  flaws.  However,  secondary  fractures,  which  were  not  observed  in  any 
of  the  experiments  on  polymers,  propagated  and  eventually  caused  the  failure  of  the  rock 
specimens.  For  example,  Fig.  2.2  shows  the  crack  growth  sequence  that  Lajtai  (1974) 
observed  in  plaster  of  paris.  What  he  labelled  as  "tensile"  fractures  appeared  first  and 
probably  correspond  to  the  wing  cracks  described  above.  Upon  further  loading,  other 
fractures  were  observed  and  were  labelled  as  "normal"  and  "inclined  shear”  fractures. 
These  labels  reflect  the  author’s  hypotheses  regarding  the  stresses  which  produced  these 
cracks. 
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Ingraffea  and  Heuze  (1980)  performed  uniaxial  compression  tests  on  limestone  and 
granodiorite  specimens  with  single  inclined  flaws.  They  reported  the  following  sequence 
of  crack  growth  (see  Fig.  2.3  for  schematic  diagram  of  observed  cracks): 

1)  Cracks  initiated  from  tensile  stress  concentration  points  near  the  flaw  tips. 
These  cracks  were  labelled  "primary"  and  correspond  to  wing  cracks  in  this  study. 

2)  These  "primary”  cracks  propagated  stably  along  a  curvilinear  path. 

3)  Another  set  of  cracks,  labelled  "secondary",  apparently  initiated  at  compressive 
stress  concentration  points  near  flaw  tips.  It  was  claimed  that  these  secondary  cracks 
actually  initiated  from  tensile  stress  concentrations  remote  from  the  flaw  tip  (see  Section 
2.2.2). 

4)  The  secondary  cracks  propagated  unstably  and  led  to  specimen  failure  at  loads 
that  were  3  to  5  times  higher  than  the  load  at  which  primary  cracks  were  first  observed. 

Wang  et  al.  (1987)  also  reported  secondary  fracturing  in  a  uniaxial  compression  test 
of  a  marble  specimen  with  a  single  inclined  flaw.  The  crack  growth  sequence  is  the  same 
as  that  observed  by  Lajtai  (1974),  and  by  Ingraffea  and  Heuze  (1980).  Tensile  cracks 
initiated  near  the  tips  of  the  flaws  and  propagated  curvilinearly  to  align  with  the  applied 
uniaxial  compressive  load.  This  was  followed  by  the  appearance  of  what  they  called  "X- 
type  shear  fractures"  which  also  intiated  from  the  flaw  tips  (see  Fig.  2.4e,  f).  These  new 
fractures  caused  the  final  rupture  of  the  specimens. 

Petit  and  Barquins  (1988)  looked  into  crack  propagation  from  single  inclined 
flaws  in  uniaxially  compressed  specimens  made  of  low  and  high  porosity  sandstones.  Fig. 
2.5  shows  a  schematic  diagram  of  the  observed  fracture  process.  They  labelled  the 
fractures  as  "bfl  for  branch  fractures,  and  "sz"  for  shear  zones.  As  in  all  the  previous 
crack  growth  sequences,  the  first  cracks  to  grow  were  the  "branch  fractures"  which  are 
equivalent  to  the  wing  cracks  in  this  study,  to  Lajtai's  and  Wang  et  al.'s  tensile  cracks,  and 
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to  Ingraffea’s  primary  cracks.  This  was  followed  by  the  appearance  of  the  "shear  zones" 
which  could  be  the  same  as  the  secondary  fracturing  that  was  observed  in  the  experiments 
described  previously. 


2.2.2  Specimens  with  Multiple  Inclined  Flaws 

In  the  early  60  s,  several  researchers  were  proposing  a  model  for  the  formation  of 
shear  faults  in  rocks  through  the  coalescence  of  pre-existing  flaws.  In  order  to  explore  the 
validity  of  this  shear  fault  model,  experiments  on  arrays  of  interacting  flaws  in  glass  and 
other  photoelastic  (polymer  )  materials  were  performed  by,  among  others,  Hoek  and 
Bieniawski  (1965),  Bombolakis  (1964)  and  Paulding  (1965).  However,  coalescence  was 
never  observed  in  these  experiments,  as  summarized  by  Brace  and  Byerlee  (1967)  in  the 
following  statement: 

"...for  each  individual  [flaw]  in  the  array,  crack  growth  followed  approximately  the 
same  directions  as  for  an  isolated  crack,  namely,  along  a  curving  path  which  approached 
the  direction  of  maximum  compression.  Clearly,  for  an  en-echelon  array,  this  would  not 
lead  to  a  through-going  fracture,  [see  Fig.  2.6]"2 1 

These  conclusions  regarding  non-coalescence  between  interacting  flaws  was 
reinforced  by  experiments  reported  by  Ingraffea  and  Heuze  (1980)  and  by  Horii  and 
Nemat-Nasser  (1986)  using  PMMA  and  CR39,  both  of  which  are  polymers.  Fig.  2.7 
shows  a  schematic  diagram  of  a  final  fracture  pattern  given  by  Ingraffea  and  Heuze  (1980). 


21  Quoted  from  W.  Brace  and  J.  Byerlee,  "Recent  experimental  studies  of  brittle  fracture  of  rocks,"  in 
Failure  and  Breakage  of  Rock,  AIME,  pp.  57-81,  1967. 
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The  figure  shows  that  wing  cracks  grew  as  in  the  single  flaw  case  but  did  not  cause  the 
coalescence  of  the  pre-existing  flaws. 

The  multiple-flaw  tests  reported  by  Horii  and  Nemat-Nasser  (1986)  were 
performed  to  establish  physical  evidence  for  a  mechanism  of  faulting  which  the  authors 
were  proposing.  Their  model  for  shear  failure  is  very  similar  to  that  being  proposed  by 
Hoek,  Brace  and  others  and  is  also  based  on  the  idea  that  a  shear  fault  is  formed  from  the 
interaction  of  many  pre-existing  flaws.  Note  that  the  original  "shear  fault"  model  stipulated 
that  the  pre-existing  flaws  would  literally  link-up  to  form  a  through-going  discontinuity. 
However,  such  coalescence  was  never  observed  in  the  previous  tests,  nor  was  it  seen  in  the 
tests  reported  by  Horii  and  Nemat-Nasser  (1986).  This  is  shown  in  Fig.  2.8  which  is  a 
schematic  diagram  of  a  plate  with  multiple,  parallel  flaws  after  it  had  been  biaxially 
compressed.  Wing  cracks  grew  out  of  the  tips  of  all  the  flaws  and  resulted  in  coalescence 
between  only  one  pair  of  pre-existing  flaws  (see  Fig.  2.8).  In  order  to  make  their  model 
for  a  shear  fault  more  plausible,  the  authors  had  to  claim  that  the  shear  fault  would  form 
due  to  the  unstable  propagation  of  wing  cracks  caused  by  interaction.  The  original  concept 
of  "linking  flaws"  had  to  be  discarded  in  view  of  the  results  of  their  model  tests. 

The  only  experiments  involving  interacting  flaws  in  compressed  rock  specimens 
were  performed  by  Wang  et  al.  (1987).  They  tested  a  marble  specimen  with  flaws  very 
similar  to  those  shown  in  Fig.  2.7.  However,  unlike  the  flaws  in  PMMA,  coalescence  did 
occur  betwen  these  flaws  in  marble  through  the  growth  of  a  crack  betwen  the  internal  flaw 
tips  (Fig.  2.9).  It  was  not  stated  whether  the  coalescing  crack  was  actually  the  intersection 
of  the  wing  cracks  growing  from  the  internal  flaw  tips,  or  whether  a  different  crack  initiated 
and  propagated  between  the  two  flaws.  The  second  possibility  seems  more  likely  in  light 
of  the  secondary  cracking  observed  for  single  flaws  in  rock  specimens  (see  Sec.  2.2.1)  as 
well  as  experiments  presented  in  this  report  (see  Section  3.3). 
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A  dependence  of  crack  growth  behavior  on  material  type  was  also  noted  in  double¬ 
edge  notch  bend  specimens  such  as  that  shown  in  Fig.  2.10  (Melin,  1989).  This  geometry 
was  being  proposed  as  a  suitable  "shear  fracture"  specimen  for  concrete2-2  since  tests 
reported  by  Bazant  and  Pffeifer  (1986)  had  shown  fracturing  in  the  region  of  high  shear 
stresses  between  the  two  notches  (see  Fig.  2.1 1).  On  the  other  hand,  tests  on  a  similar 
specimen  geometry  but  using  PMMA  show  cracks  that  initiate  at  the  edge  of  the  notches 
and  propagate  in  a  curvilinear  path  (see  Fig.  2.12).  The  difference  between  the  fracture 
behavior  in  Figs.  2.1 1  and  2.12  is  similar  to  the  observed  coalescence  and  non-coalescence 
between  flaws  in  rock  and  polymers. 


2  2  Concrete  and  rock  arc  both  classified  as  cohesive,  frictional  materials.  Thus,  the  behavior  of  these 
materials  is  expected  to  be  similar. 
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Figure  2. 1  Typical  crack  growth  from  a  pre-existing  flaw  in  polymer  specimens 
subjected  to  uniaxial  compression. 
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PRE-EXISTING  ELASTIC  FLAW  -  TENSILE  FRACTURE 
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Figure  2.3 


Crack  growth  from  a  single  inclined  flaw  in  a  limestone  specimen  under 
uniaxial  compression  [from  Ingraffea  and  Heuze,  1980]. 
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high  porosity  sandstone 


low  porosity  sandstone 


Figure  2.5  Crack  growth  from  a  single-inclined  flaw  in  low  and  high  porosity 

sandstone  specimens  under  uniaxial  compression  [from  Petit  and  Barquins, 
1988]. 
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Figure  2.6 


Crack  growth  from  multiple,  en-echelon  flaws  in  a  glass  specimen  under 
uniaxial  compression  [from  Brace  and  Byerlee,  1967], 
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Figure  2.7 


Figure  2.8 


Crack  growth  from  multiple  flaws  in  a  PMMA  specimen  under  uniaxial 
compression  [from  Ingraffea  and  Heuze,  1980]. 


Crack  growth  from  multiple  flaws  in  a  CR39  specimen  under  uniaxial 
compression  [from  Horii  and  Nemat-Nasser,  1986]. 
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2.3  Modelling  of  Crack  Growth  from  Inclined  Flaws 


There  are  two  general  approaches  to  modelling  crack  growth:  (1)  the  discrete  crack 
approach  and  (2)  the  smeared  crack  approach.  In  the  discrete  crack  approach  (e.g. 
Ingraffea  and  Heuze,  1980;  Chan,  1986;  Annigeri,  1984),  certain  rules  are  applied  to 
determine  the  direction  and  load  at  which  crack  extension  occurs.  These  rules  are  called 
crack  propagation  criteria,  examples  of  which  are  given  in  Table  2.1. 


Table  2.1  Summary'  of  Crack  Propagation  Criteria 

Erdogan  and  Sih  (1963) 

Maximum  Circumferential 

Stress  Criterion 

Crack  propagates 
perpendicular  to  the 

direction  of  maximum 

tensile  hoop  stress. 

Kipp  and  Sih  (1974) 

Minimum  Strain  Energy 
Density  Criterion 

Crack  propagates  along  the 

direction  of  minimum  strain 

energy  density 

Chang (1974) 

Maximum  Circumferential 

Strain  Criterion 

Crack  propagates 
perpendicular  to  the 
direction  of  maximum 

tensile  hoop  strain. 

Once  the  direction  has  been  established  using  any  one  of  these  criteria,  an  actual 
discontinuity  is  attached  to  the  existing  crack.  Subsequent  crack  extension  is  determined 
from  the  stress/strain  energy/strain  field  at  the  tip  of  this  added  discontinuity. 
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In  the  smeared  crack  approach  (e.g.  Bazant  and  Oh,  1983;  Rots  et  al.,  1985), 
fracturing  is  modelled  through  a  constitutive  law  that  simulates  the  effects  of  material 
degradation  (e.g.  loss  of  stiffness  and  strain  softening).  Thus,  if  the  constitutive  law 
assumes  that  material  degradation  is  dependent  on  tensile  principal  stresses,  then  the 
"fractured”  zone  would  be  expected  to  lie  along  the  most  highly  stressed  regions  but  would 
not  necessarily  be  perpendicular  to  the  tensile  principal  stresses  as  in  the  discrete  crack 
approach. 

Both  the  discrete  crack  and  the  smeared  crack  approaches  have  been  applied  to 
modelling  mixed-mode  crack  growth  in  notched  three-point  and  four-point  bend  specimens 
(Ingraffea  and  Panthaki,  1985;  Swartz  and  Taha,  1990;  Rots,  1986).  However,  only  the 
discrete  crack  approach  has  been  used  to  simulate  crack  growth  from  an  inclined  flaw  in 
uniaxial  compression  (Ingraffea  and  Heuze,  1980). 

Ingraffea  and  Heuze  (1980)  employed  the  finite  element  method  and  applied 
existing  crack  propagation  criteria  to  successfully  simulate  the  "primary"  (wing)  crack 
growth  which  they  had  observed  in  their  experiments.  The  criteria  which  they  used  were: 
the  maximum  circumferential  stress  criterion  (Erdogan  and  Sih,  1963)  and  the  minimum 
strain  energy  density  criterion  (Kipp  and  Sih,  1975),  see  Table  2.1.  A  comparison  of  the 
wing  crack  paths  are  shown  in  Fig.  2.13.  Note  that  the  simulated  paths  follow  the  general 
shape  of  the  observed  primary  (wing)  crack,  i.e.,  the  simulated  paths  are  curvilinear  and 
eventually  align  with  the  applied  compressive  load.  In  explaining  the  growth  of  the 
secondary  cracks,  Ingraffea  and  Heuze  looked  at  the  stress  field  after  the  growth  of  the 
wing  crack  and  found  a  region  of  tensile  principal  stresses  around  the  location  of  the 
secondary  crack  (see  Fig.  2.14).  Since  rock  is  weak  in  tension,  they  claimed  that  the 
secondary  crack  initiated  within  this  tensile  region,  specifically  at  point  A  (Fig.  2.14)  which 
is  at  a  distance  from  the  notch  tip.  Once  initiated,  secondary  crack  growth  was  driven  by 
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the  tensile  principal  stresses.  The  possibility  of  the  secondary  crack  initiating  from  the 
point  of  compressive  stress  concentration  at  the  notch  tip  was  dismissed  on  the  basis  of 
photographic  proof  which,  however,  was  not  explicitly  described. 

Only  Horii  and  Nemat-Nasser  (1986)  have  attempted  to  analytically  model  crack 
propagation  from  multiple  flaws  in  compression.  Wing  cracks  were  represented  by  a 
superposition  of  dislocations  such  that  when  boundary  conditions  were  enforced, 
dislocation  densities  could  be  determined  from  a  system  of  integral  equations  (for  more 
details,  refer  to  Horii  and  Nemat-Nasser,  1986).  Interaction  between  adjacent  flaws  was 
approximated  by  the  "method  of  pseudotractions"  which  Horii  and  Nemat-Nasser  (1985) 
had  also  developed.  This  method,  however,  can  only  be  used  for  an  infinite  line  of  parallel 
and  equidistant  flaws  (see  Fig.  2.15).  Furthermore,  since  their  experiments  in  CR39  only 
showed  this  kind  of  fracturing,  wing  crack  growth  alone  was  modelled  and  coalescence 
was  not  observed  in  their  simulations. 
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Figure  2. 1 3  Comparison  of  observed  and  simulated  primary  cracks,  simulated  crack  (1 ) 

generated  using  Maximum  Circumferential  Stress  Criterion,  simulated  crack 
(2)  generated  by  Minimum  Strain  Energy  Density  Criterion  [from  Ingraffea 
and  Heuze,  1980], 


Figure  2.14  Tensile  stress  region  (marked  by  arrows)  near  tips  of  a  pre-existing  flaw 
with  a  primary  crack.  Secondary  crack  (see  Fig.  2.13)  thought  to  have 
initiated  at  pt.  A  [from  Ingraffea  and  Heuze,  1980]. 
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Figure  2.15  Shear  fault  formation  through  crack  propagation  from  multiple  pre-existing 
flaws.  P-P'  and  similar  segments  are  pre-existing  flaws.  P'-Q'  and  P-Q 
are  wing  cracks  [from  Horii  and  Nemat-Nasser,  1980]. 


Chapter  3 

Experimental  Procedures  and  Results 


3.1  Introductory  Remarks 

As  pointed  out  in  Chapter  2,  crack  propagation  criteria  from  fracture  mechanics 
have  been  successfully  used  to  simulate  wing  crack  growth  observed  in  compressed 
polymer  blocks  with  pre-existing  flaws  (Horii  and  Nemat-Nasser,  1985;  Ingraffea  and 
Heuze,  1980).  On  the  other  hand,  initial  experiments  on  flawed  gypsum  blocks  exhibited 
fracture  behavior  that  was  quite  different  from  what  was  observed  in  polymer  materials  and 
predicted  by  fracture  mechanics  theory.  In  particular,  coalescence  occured  through  the 
growth  of  secondary  cracks  that  were  distinct  from  the  wing  cracks. 

There  was  thus  a  discrepancy  between  different  experimental  and  analytical  studies, 
and  a  systematic  investigation  was  necessary  to  determine  coalescence  mechanisms  of 
flaws  in  rock-like  materials.  This  chapter  presents  the  experimental  aspect  of  this 
investigation.  It  is  subdivided  into  two  main  sections:  Section  3.2  contains  the  procedures 
for  specimen  fabrication,  loading  and  crack  monitoring  while  Section  3.3  concentrates  on 
the  results  of  the  experiments  in  form  of  photographs  and  loads  corresponding  to  the 
various  stages  of  the  crack  growth  process.  Some  observations  relating  to  the  surface 
characteristics  of  the  wing  and  secondary  cracks  are  also  included  in  Section  3.3. 
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3.2  Experimental  Procedures 


3.2.1  Model  Material 


The  material  chosen  for  this  study  is  hydrated  gypsum  made  from  a  mixture  of 
Hydrocal  B-l  1  (manufactured  by  US  Gypsum  Company),  celite  and  water.  The  following 
mass  proportions  were  used: 


water 

gypsum 


0.40 


water  _  .. 
celite  "  " 


Approximate  tensile  and  compressive  strengths  were  obtained  from  data  collected 
by  Nelson  (1968)  on  a  mix  that  is  similar  to  what  was  used  in  this  study: 


water 

gypsum 


0.45 


water  _ 
celite 


Uniaxial  Compressive  Strength  =  23  MPa 


Uniaxial  Tensile  Strength  =  2.3  MPa. 

Thus,  this  mix  exhibits  a  compressive  strength/tensile  strength  ratio  of  about  10.0,  which  is 
in  the  lower  range  of  ratios  for  naturally  occuring  rocks. 

Gypsum  was  selected  as  a  model  material  primarily  due  to  the  ease  with  which 
flaws  can  be  introduced  into  the  specimens  by  pre-moulding  them  in.  Other  natural  rocks 
were  considered,  but  the  costs  for  cutting  harder  rocks  to  simulate  flaws  would  have 
restricted  the  number  of  flaw  arrangements  tested. 
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3.2.2  Specimen  Geometry 


The  dimensions  of  the  gypsum  blocks  are  shown  in  Figure  3.1.  The  thickness  of 
the  blocks  varied  between  28  to  30  mm.  The  specimen  length  was  set  at  twice  the  width  in 
order  to  reduce  end  effects.  Parallel  flaws  that  were  12.7  mm  (0.5")  long  and 
approximately  0.254  mm  (0.01")  wide  were  molded  into  the  blocks  such  that  the  planes  of 
the  flaws  were  perpendicular  to  the  76.2  mm  x  152.4  mm  face  of  the  the  block.  The  model 
tests  were  thus  limited  to  the  investigation  of  the  two-dimensional  fracture  patterns 
associated  with  these  flaws.  It  was  also  assumed  that  the  out-of-plane-stress  had  a  minimal 
effect  on  the  fracturing. 

The  geometric  parameters  which  characterize  the  flaw  arrangements,  namely  the 
flaw  inclination  angle  and  the  ligament  inclination  angle,  are  shown  in  Fig.  3.2.  The 
"external"  and  "internal”  flaw  tips  are  also  identified  in  this  figure  since  these  terms  will  be 
used  in  the  subsequent  descriptions  of  crack  growth  (Section  3.3).  The  flaw  arrangements 
could  be  classified  according  to  non-overlapping  and  overlapping  flaws  (Fig.  3.3),  where 
the  non-overlapping  flaws  are  characterized  by  a  ligament  inclination  angle  which  is  less 
than  or  equal  to  90',  while  the  overlapping  flaws  have  ligament  inclination  angles  that  are 
greater  than  90'.  The  overlap  is  determined  with  respect  to  the  direction  of  the  applied 
uniaxial  compression  which  coincides  with  the  longer  axis  of  the  block  (Fig.  3.3).  Three 
inclination  angles  were  used  in  the  experiments:  30',  45'  and  60'.  The  non-overlapping 
geometries  for  each  inclination  angle  are  shown  in  Figs.  3.4  through  3.6  while  the 
overlapping  geometries  for  all  three  inclination  angles  are  shown  in  Fig.  3.7.  In  the 
following  discussion,  the  following  convention  for  labelling  flaw  arrangments  will  be  used: 
"X'-Y'M  where  X  is  the  flaw  inclination  angle  and  Y  is  the  ligament  inclination  angle.  For 
example,  ”30'-45'"  will  refer  to  flaws  inclined  at  30'  having  a  ligament  that  is  inclined  at 
45'  with  respect  to  the  horizontal. 
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All  flaw  arrangements  had  ligament  lengt;  s  of  12.7  mm  which  is  equal  to  the  length 


of  the  flaws. 


3.2.3  Specimen  Fabrication 

The  mold  consisted  of  4  stainless  steel  plates  that  were  bolted  together  to  form  a 
76.2mm  x  152.4  mm  block.  These  four  plates  were  fastened  onto  a  stainless  steel  base 
plate  that  was  then  clamped  to  a  vibrating  table.  The  inside  surfaces  of  the  four  steel  plates 
were  wrapped  in  vinyl  sheets  to  prevent  bubbling  at  the  steel  shims  which  was  observed 
when  the  mold  was  left  unlined.31 

The  flaws  were  pre-molded  by  steel  shims  which  were  arranged  in  the  mold  to 
produce  the  desired  flaw  geometries.  The  shims  were  pulled  out  after  the  gypsum  had 
hardened  sufficiently.  The  step-by-step  procedure  for  mold  preparation  was  as  follows 
(see  Figure  3.8): 

(1)  A  stiff  piece  of  cardboard  was  cut  to  the  dimensions  of  the  block  (76.2mm  x  152.4 
mm).  This  cardboard  was  wrapped  in  a  vinyl  sheet  to  prevent  it  from  absorbing 
water  from  the  pre-hardened  gypsum.  Two  slits  were  cut  into  the  cardboard, 
corresponding  to  the  desired  flaw  arrangements. 


31  The  bubbling  could  be  related  to  an  "electrolytic”  current  between  the  stainless  steel  plates  and  the  non- 
stainlcss  steel  shims.  It  was  oberved  that  the  bubbling  was  significantly  reduced  when  the  steel  plates 
where  wrapped  in  self-adhesive  vinyl  shelf  paper. 
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(2)  Filler  gage  steel  stock  with  a  width  of  12.7  mm  (0.5”)  and  thickness  of  0.10  mm 
(0.004")  was  cut  into  50.8  mm  (2")  long  pieces.  Two  of  these  steel  shims  were 
inserted  into  qq£  of  the  slits  in  the  cardboard.  Two  0.004"  shims  were  used  rather 
than  one  0.008"  shim  because  it  was  more  difficult  to  pull  out  one  0.008"  shim  than 
it  was  to  pull  out  one  0.004"  shim  at  a  time.  The  shims  were  also  greased  to 
prevent  the  gypsum  from  adhering  to  the  steel. 

(3)  Two  slits  were  cut  into  a  shorter  piece  of  76.2  mm  wide  cardboard.  This  smaller 
piece  was  used  to  ensure  that  the  shims  were  more  or  less  perfectly  vertical  while 
the  gypsum  was  being  poured  into  the  mold.  (See  Figure  3.8.) 

(4)  The  entire  cardboard-steel  shim  assembly  was  then  set  into  the  stainless  steel  mold. 

In  order  to  have  blocks  with  uniform  properties,  the  mixing  and  curing  process  was 

strictly  controlled.  The  mixing/curing/polishing  procedure  was  as  follows: 

( 1 )  The  water  and  celite  were  combined  and  mixed  for  30  seconds. 

(2)  The  gypsum  was  gradually  added  to  the  water-celite  mixture  and  blended  for  four 
minutes. 

(3)  The  mixture  was  immediately  poured  into  the  mold  after  blending.  The  mold  was 
vibrated  during  and  a  few  minutes  after  the  mixture  was  poured  into  it  in  order  to 
remove  entrained  bubbles. 

(4)  When  the  mixture  had  sufficiently  hardened,  the  blocks  were  taken  out  of  the  mold. 
The  shims  were  pulled  out  prior  to  polishing  and  curing. 
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(5)  The  blocks  were  allowed  to  air  dry  for  about  a  day,  after  which  the  76.2  x  152.4 
mm  surfaces  were  ground  and  polished.  This  made  it  easier  for  fractures  to  be 
observed  during  testing. 

(6)  The  mass  of  the  blocks  was  measured  before  they  were  cured  in  an  oven  at  40 °C. 
The  blocks  were  occasionally  weighed  during  the  curing  process.  When  the  mass 
had  stabilized  indicating  that  all  free  water  had  been  removed,  the  blocks  were  ready 
for  testing. 
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3.2.4  Loading  and  Crack  Monitoring  Procedure 


Specimens  were  uniaxially  compressed  in  a  computer-controlled  INSTRON  Model 
1331  servohydraulic  machine.  The  compression  loading  fixture  had  a  ball  and  socket  joint 
that  compensated  for  specimen  ends  that  were  not  perfecdy  parallel.  Step  loading  was  done 
under  displacement  control  with  -0.05  -  0.15  mm  per  step3-2  applied  at  a  rate  of  0.005 
mm/sec.  After  each  load  step  was  applied  and  without  unloading  the  specimen,  the 
specimen  surfaces  were  scanned  for  newly  formed  fractures  through  a  low  powered 
microscope  (20x).  The  stereomicroscope  Lad  to  be  equipped  with  a  stand  that  allowed  one 
to  view  the  specimen  surface  while  translating  along  a  vertical  plane. 

Early  tests  were  performed  with  continuous  loading  and  a  video  camera  to  record 
crack  growth  [Reyes  et  al.  1989];  however,  the  video  equipment  used  did  not  have  enough 
resolution  to  pick  up  some  of  the  hair-line  cracks.  The  crack  growth  at  failure  was  also 
very  rapid  when  the  specimens  were  continuously  loaded.  It  was  found  that  some  cracks, 
which  were  unobservable  with  the  naked  eye,  were  clearly  visible  with  the  microscope. 
Furthermore,  step  loading  was  adopted  since  compression  crack  growth  is  usually  stable 
prior  to  final  specimen  failure  and  this  allowed  more  time  to  scan  the  critical  regions  on  the 
specimen  surface.  Step  loading  and  crack  growth  monitoring  through  a  microscope 
resulted  in  a  better  documentation  of  the  crack  growth  sequence  during  uniaxial 
compression. 


3  2  Lower  range  is  used  when  block  is  close  to  failure  point. 
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Figure  3.1  Geometry  of  prismatic  specimens  with  pre-existing  flaws. 
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Figure  3.2  Parameters  defining  flaw  geometry. 
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Figure  3.3  Non-overlapping  (left)  and  overlapping  (right)  flaws. 
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L.I.A.  =  45*  L.I.A.  =  60'  L.I.A.  =75’  L.I.A 


Mold  preparation  procedure. 


3.3  Experimental  Results 


3.3.1  Sequence  and  Patterns  of  Fracture  Growth 

For  all  flaw  geometries,  wing  cracks  at  the  external  and  internal  flaw  tips  appeared 
first.  Subsequent  fracture  behavior  was  more  dependent  on  the  flaw  geometry,  particularly 
on  whether  the  flaws  were  overlapping  or  not: 

(1 )  For  all  non-overlapping  flaw  geometries,  the  internal  wing  crack  growth  was  arrested 
(i.e.,  wing  cracks  did  not  propagate  through  the  entire  specimen)  and,  with  an  increase  in 
load,  was  followed  by  the  appearance  of  secondary  cracks.  These  secondary  cracks  caused 
coalescence  between  flaws  in  all  non-overlapping  arrangements  except  between  45'-45* 
and  30°-30°  flaws  (note  that  both  are  coplanar  geometries). 

(2)  For  the  overlapping  flaw  geometries,  coalescence  occured  through  wing  crack  growth 
for  most  cases.  There  was  one  case  where  it  occured  through  secondary  fracturing  and  two 
cases  where  coalescence  did  not  take  place  at  all  prior  to  complete  disintegration  of  the 
specimen. 


Coalescence  between  Non-overlapping  Flaws 

Fracturing  between  non-overlapping  flaws  was  very  consistent  in  specimens  with 
identical  flaw  geometries,  i.e.,  tests  were  reproducible.  Furthermore,  except  for  two 
coplanar  cases  which  will  be  described  later,  secondary  crack  growth  resulted  in  the 
coalescence  of  the  internal  flaw  tips.  Differences  among  coalescence  mechanisms  for 
various  non-overlapping  flaw  geometries  were  found  in:  (1)  measured  loads  at  which 
coalescence  occured  (Section  3.3.2)  and  (2)  surface  characteristics  of  coalescence  cracks 
(Section  3.3.3).  However,  since  the  sequence  of  crack  propagation  was  generally  the  same 
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for  all  the  non-overlapping  geometries,  only  the  crack  growth  processes  for  three 
representative  flaw  arrangements  will  be  described,  namely  30*-45\  45'-90°  and  60°-60° 
flaws. 


Crack  growth  from  30  - 45  ° flaws 

Fig.  3.9  is  a  load-displacement  curve  obtained  from  uniaxial  compression  of  a 
block  with  30°-45°  flaws.  The  loads  at  which  the  photographs  in  Figs.  3.10  through  3.12 
were  taken  are  also  indicated  on  this  graph.  Figure  3.10(a)  and  (b),  taken  at  an  applied 
load  of  47  kN  (21  MPa)  3  3  ,  show  wing  cracks  propagating  out  of  the  flaw  tips.  The 
wing  cracks  were  first  observed  at  a  load  of  25  kN  (1 1.2  MPa).  Slight  spalling  or  crushing 
is  also  evident  in  both  Figs.  3.10(a)  and  (b)  as  marked  by  the  arrows.  At  an  applied  load 
of  54  kN  (24.1  kN),  coalescence  suddenly  occured  through  the  growth  of  a  secondary 
crack  that  connected  the  internal  flaw  tips  (see  Fig.  3.11).  Note  that  the  previously 
observed  crushing/spalling  became  a  part  of  this  coalescence  crack;  however,  the 
coalescence  crack  propagated  so  rapidly  that  it  was  impossible  to  establish  whether  it 
propagated  from  the  center  of  the  ligament  or  from  the  flaw  tips.  Fig.  3.12  is  an  expanded 
view  of  the  flaws  after  coalescence.  The  external  wing  cracks  had  opened  and  extended  at 
this  point.  Whether  the  external  wing  cracks  extended  before  or  after  coalescence  could  not 
be  established  due  to  the  rapidity  of  crack  growth.  However,  it  seems  likely  that  the 
coalescence  crack  preceded  the  opening  and  extension  of  the  external  wing  cracks.  What 


3  3  In  the  subsequent  discussion,  the  applied  force  will  be  given  followed  by  the  equivalent  applied  traction 
enclosed  in  parentheses.  The  applied  traction  is  equal  to  the  applied  force  divided  by  the  specimen's  loading 
surface. 
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could  have  happened  was  the  following:  (1)  the  linkage  of  the  flaws  resulted  in  a  longer 
"effective  flaw",  (2)  the  increase  in  "effective  flaw"  length  increased  the  relative  sliding 
displacements  of  the  flaw  surfaces,  (3)  the  increase  in  relative  sliding  displacements  caused 
the  external  wing  cracks  to  open,  (4)  the  opening  of  the  wing  cracks  then  increased  the 
singularity  at  the  wing  crack  tips  and  drove  the  wing  cracks  to  lengthen. 
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Load  (kN) 


Displacement  (mm) 


Figure  3.9  Load-displacement  curve  taken  from  uniaxial  compression  of  a  gypsum 

block  with  30°-45°  flaws.  Loads  at  which  subsequent  photos  were  taken  are 
marked  on  the  graph. 
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Figure  3.10  (a)  and  (b)  Wing  crack  growth  and  spalling  from  internal  30* -45*  flaw  tips 

[applied  load  =  47kN]. 
sp  =  spalling;  int.  wc  =  internal  wing  crack. 


Figure  3.11 


Secondary  crack  growth  resulting  in  coalescence  between  internal  30*-45* 
flaw  tips  [applied  load  =  54kN], 
sc  =  secondary  crack. 
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Figure  3. 12  Expanded  view  of  30’-45*  flaws  after  coalescence.  Note  open  and  extended 
external  wing  cracks  [applied  lead  =  54kN]. 
sc  =  secondary  cracks;  ext  wc  =  external  wing  crack. 
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Crack  growth  from  45  °-90  ° flaws 


Figure  3.13  is  a  load-displacement  curve  obtained  from  uniaxial  compression  of  a 
block  with  45*-90*  flaws.  The  loads  at  which  the  photographs  in  Figures  3.14  through 
3.18  were  taken  are  also  indicated  on  this  graph.  In  this  test,  the  wing  cracks  were  first 
observed  at  an  applied  load  of  27  kN  (12.5  MPa).  Wing  crack  growth  with  slight  spalling 
at  the  internal  flaw  tips  are  shown  in  Figs.  3.14  taken  at  applied  loads  of  37  kN  (17.1 
MPa).  At  a  load  of  50  kN  (23.2  MPa),  the  external  wing  cracks  [Figs.  3.15(a)  and  (c)] 
were  longer  than  the  internal  wing  cracks  [Figs.  3.15(b)],  i.e.,  the  internal  wing  cracks 
were  arrested  upon  reaching  the  level  of  the  adjacent  flaw  tip.  At  a  load  of  51  kN  (26.4 
MPa),  spalling  increased  at  the  internal  flaw  tips  (see  Fig.  3.16)  and  finally,  at  a  load  of  61 
kN  (28.3  MPa),  coalescence  occured  through  the  sudden  growth  of  a  crack  in  the  zone 
between  the  internal  wing  cracks  (see  Fig.  3.17).  An  expanded  view  of  the  final  fracture 
pattern  after  coalescence  is  shown  in  Fig.  3.18.  As  in  the  30*-45*  case,  the  external  wing 
cracks  opened  and  extended  while  the  internal  wing  cracks  closed  after  coalescence. 
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Figure  3. 1 3  Load-displacement  curve  taken  from  uniaxial  compression  of  a  gypsum 

block  with  45°-90°  flaws.  Loads  at  which  subsequent  photos  were  taken  are 
marked  on  the  graph. 
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Figure  3.14  Wing  crack  and  s palling  from  internal  45*-90*  flaw  tips  [applied  load  = 
37kN] 

sp  =  spalling;  int.  wc  =  internal  wing  cracks. 
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int.  wc 
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Figure  3.17  Secondary  crack  growth  resulting  in  coalescence  between  internal  45*-90* 
flaw  tips  [applied  load  =  61kN]. 
sc  =  secondary  crack 
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Figure  3.18  Expanded  view  of  45*-90*  flaws  after  coalescence.  Note  open  and  extended 
external  wing  cracks  [applied  load  =  61kN]. 
sc  =  secondary  crack,  ext  wc  =  external  wing  crack 
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Crack  growth  from  60  ° -60  ° flaws 


Figure  3.19  is  a  load-displacement  curve  obtained  from  uniaxial  compression  of  a 
block  with  60*-60‘  flaws.  The  loads  at  which  the  photographs  in  Figure  3.20  through 
3.223-4  were  taken  are  also  indicated  on  this  graph.  Wing  cracks  were  first  visible  at  an 
applied  load  of  56  kN  (25.3  MPa).  At  70  kN  (31.7  MPa),  spalling  or  possibly  secondary 
cracking  initiated  from  the  internal  tips  (Fig.  3.20)  but  did  not  "connect"  until  the  applied 
load  was  increased  to  73  kN  (33  MPa),  see  Fig.  3.21.  This  test  was  one  of  the  few  cases 
where  coalescence  clearly  intiated  from  the  flaw  tips.  An  expanded  view  of  the  flaws  is 
shown  in  Fig.  3.22.  Though  not  visible  in  the  photograph,  the  external  wing  cracks  had 
opened  and  extended  after  coalescence. 


3  4  In  this  test,  red  ink  was  used  to  mark  the  observed  cracks. 
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Displacement  (mm) 


Figure  3. 19  Load-displacement  curve  taken  from  uniaxial  compression  of  a  gypsum 

block  with  60*-60*  flaws.  Loads  at  which  subsequent  photos  were  taken 
are  marked  on  the  graph. 
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Figure  3.21  Coalescence  between  internal  60*-60*  flaw  tips  [applied  load  =  73kN]. 
sc  =  secondary  crack 
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Figure  3.22  Expanded  view  of  60’-60*  flaws  after  coalescence,  external  wing  cracks 
opened  and  extended  at  this  load  [applied  load  =  73kN]. 
sc  =  secondary  crack,  ext  wc  =  external  wing  crack 
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Crack  growth  from  flaws  with  other  non-overlapping  geometries 


Coalescence  crack  patterns  for  other  non-overlapping  flaw  geometries  are  shown  in 
Figs.  3.23  through  3.28.  Fig.  3.23  shows  the  coalescence  crack  for  a  specimen  with  30*- 
45*  flaws,  identical  to  the  one  shown  in  Fig.  3.11  (also  shows  30*-45*  flaws).  Figs.  3.24 
and  3.25  are  photographs  of  coalescence  cracks  for  two  identical  specimens  with  30*-75‘ 
cracks.  The  similarity  between  the  coalescence  cracks  in  Figs.  3.24  and  3.25  is  rather 
striking;  both  are  oriented  parallel  to  the  applied  vertical  compressive  load.  Figs.  3.26  and 
3.27  show  coalescence  cracks  for  identical  specimens  with  45*-60*  flaws.  One  can 
observe  that  the  coalescence  cracks  between  45*-60‘  flaws  consist  of  "shearing"  and 
"tensile"  segments,  as  indicated  in  the  figures.  Fig.  3.28  shows  the  coalescence  crack 
between  60°-45*  flaws.  It  is  interesting  to  see  that  coalescence  took  place  for  this  geometry 
even  though  the  flaws  were  left-stepping  (see  schematic  drawing  on  figure)  rather  than 
right  stepping  like  the  other  geometries. 

In  most  of  the  specimens  tested,  the  external  wing  cracks  did  not  propagate  all  the 
way  to  the  top  and  bottom  edges  of  the  specimens.  These  specimens  remained  intact  after 
they  were  unloaded,  and  in  order  to  inspect  the  surface  of  the  coalescence  cracks  (see  Sec. 
3.3.3  for  a  discussion  of  surface  characteristics),  they  had  to  be  cut  as  shown  in  Fig.  3.29. 

This  process  of  cutting  was  useful  for  verifying  or  refuting  observed  coalescence 
between  pre-existing  flaws.  For  example,  Fig.  3.30  shows  an  "apparent"  coalescence 
crack  that  connected  the  internal  tips  of  two  45*-45*  flaws  at  an  applied  load  of  64.4  kN  (29 
MPa) .  When  this  crack  appeared,  the  specimen  was  unloaded  because  it  was  thought  that 
linkage  had  already  taken  place.  However,  upon  cutting  this  specimen,  it  could  not  be 
pulled  apart  along  the  ligament  This  clearly  indicated  that  coalescence  had  not  occured  and 
that  the  "apparent"  coalescence  crack  was  actually  an  initial  stage  of  spalling  within  the 
ligament  area. 
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Fig.  3.30  also  shows  that  the  external  wing  cracks  remained  closed  even  after  the 
appearance  of  the  "apparent"  coalescence  crack.  Recalling  the  earlier  explanation  regarding 
the  opening  of  external  wing  cracks  as  an  effect  of  coalescence,  this  observation  is  not  at  all 
surprising  since  the  45'-45*  flaws  were  not  connected  at  the  time  the  photo  was  taken.  It 
was  also  observed  from  the  load-displacement  curves  in  Fig.  3.9,  3.13  and  3.19  that, 
under  displacement  control,  coalescence  was  always  accompanied  by  a  drop  in  the  load. 
Such  a  load  drop  was  not  present  in  the  load-displacement  curve  of  the  specimen  with  45'- 
45'  flaws.  This  is  also  consistent  with  the  fact  that  coalescence  did  not  take  place  between 
the  flaws  in  this  specimen. 

The  "real"  coalescence  behavior  of  45'-45*  flaws  could  not  be  obtained  from  the 
specimen  shown  in  Fig.  3.30.  Unfortunately,  the  second  specimen  with  45'-45*  flaws 
disintegrated  completely  when  the  testing  machine  malfunctioned  and  overloaded  the  block. 
Hence,  the  fracture  coalescence  behavior  for  45*-45‘  flaws  could  not  be  established  by  the 
test  series  performed  in  this  study. 

The  fracture  behavior  of  30'-30‘  flaws  was  also  different  from  the  rest  of  the  non¬ 
overlapping  flaw  geometries.  A  specimen  with  30'-30’  flaws  is  shown  in  Fig.  3.31  at  an 
applied  load  of  64  kN  (29  MPa).  Substantial  spalling  can  be  seen  at  the  flaw  tips,  and 
when  failure  occured  at  67  kN  (30  MPa),  the  secondary  cracks  did  not  connect  the  flaw 
tips  (see  Fig.  3.32)  as  they  did  in  the  other  flaw  geometries.  The  second  specimen  with 
30*-30*  flaws  was  loaded  until  spalling  similar  to  that  shown  in  Fig.  3.31  occured.  The 
specimen  was  then  unloaded,  and  an  attempt  to  section  the  ligament  area  was  made.  It  was 
thought  that  sections  taken  immediately  prior  to  failure  would  show  some  microstructural 
evidence  of  progressive  fracturing.  However,  gypsum  has  a  strong  tendency  to  pulverize 
and  suffered  too  much  disturbance  after  the  sectioning  process.  Thus,  no  evidence  of  pre¬ 
fracturing  could  be  observed. 
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Figure  3.23  Secondary  crack  growth  resulting  in  coalescence  between  30*-45  flaws 
compare  with  Fig.  3.11. 
sc  =  secondary  crack 


Figure  3.24  Secondary  crack  growth  resulting  in  coalescence  between  30* -75*  flaws, 
compare  with  Fig.  3.25. 
sc  =  secondary  crack 
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Figure  3.25  Secondary  crack  growth  resulting  in  coalescence  between. 30* -75*  flaws, 
compare  with  Fig.  3.24. 
sc  =  secondary  crack 
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Figure  3.26  Secondary  crack  growth  resulting  in  coalescence  between  45*-60*  flaws, 
specimen  #1  (compare  with  Fig.  3.27). 
sc  =  secondary  crack 


Figure  3.27  Secondary  crack  growth  resulting  in  coalescence  between  45*-60*  flaws, 
specimen  #2  (compare  with  Figure  3.26). 
sc  =  secondary  crack 
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Figure  2 


Figure  3.30  "Apparent"  coalescence  crack  (marked  by  arrows)  between  45**45'  flaws. 
£  Closed  external  wing  cracks  hardly  visible  in  photo. 
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Figure  3.3 1  Spalling  around  internal  30*-30*  flaw  tips  [applied  load  =  64kN]. 
sp  =  spalling 
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Figure  3.32  Secondary  crack  growth  without  coalescence  between  SO^O*  flaws 
[applied  load  =  67kN]. 

int  wc  =  internal  wing  crack,  ext.  wc  =  external  wing  crack,  sc  = 
secondary  crack. 


Fig.  3.33(a)  and  (b)  are  photographs  from  separate  tests  on  specimens  with  30*- 
120*  flaws  under  applied  loads  of  62.4  kN  (27.5  MPa)  and  61  kN  (27.4  MPa), 
respectively.  Linkage  of  the  flaws  is  not  evident  in  the  photographs,  but  cutting  both 
specimens  showed  that  in  each  case  the  internal  wing  cracks  were  connected  to  the  adjacent 
internal  flaw  tips,  i.e.,  coalescence  occured  through  wing  crack  growth.  Coalescence  was 
also  indicated  by  drops  in  load  under  displacement  control  and  by  the  open  external  wing 
cracks  (see  discussion  on  coalescence  in  non-overlapping  flaws).  The  fracture  behavior 
was  also  remarkably  consistent  between  the  two  identical  specimens. 

For  one  of  the  specimens  with  30*- 135*  flaws  (Fig.  3.34),  coalescence  occured 
through  the  growth  of  the  lower  internal  wing  crack  towards  the  upper  internal  flaw  tip  at 
an  applied  load  of  67.5  kN  (29.2  MPa).  Coalescence  is  rather  obvious  from  the 
photograph,  and  is  also  indicated  by  the  open  external  wing  cracks.  Significant  spalling 
and  other  cracks  (marked  by  "A"  in  Fig.  3.34)  occured  simultaneously  with  coalescence. 
The  other  specimen  with  30°-135*  flaws  is  shown  in  Fig.  3.35  under  an  applied  load  of 
53.5  kN  (24.7  MPa).  The  wing  cracks  from  both  internal  flaw  tips  appear  to  have 
interesected  to  form  a  coalescence  crack.  However,  this  "apparent"  coalescence  crack  as 
well  as  the  external  wing  cracks  were  closed  under  the  applied  load  of  53.5  kN.  Loading 
of  this  specimen  was  discontinued  shortly  after  the  photograph  in  Fig.  3.35  was  taken 
because  of  splitting  along  the  external  edge  of  the  block.  Upon  cutting  (as  in  Fig.  3.29), 
the  block  had  to  be  pulled  apart  with  some  force  along  the  "apparent"  coalescence  crack, 
thus  establishing  that  the  crack  shown  in  Fig.  3.35  did  not  actually  connect  the  internal  flaw 
tips. 


In  one  of  the  specimens  with  45*- 135*  flaws  (see  Fig.  3.36)  under  an  applied  load 
of  72  kN  (33  MPa),  secondary  crack  growth  connected  each  external  flaw  tip  with  the 
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internal  tip  of  the  adjacent  flaw.  This  was  accompanied  by  the  opening  of  the  external  wing 
cracks  and  by  a  load  drop  under  displacement  control.  In  the  other  specimen,  the  lower 
internal  wing  crack  suddenly  extended  toward  the  upper  flaw  at  an  applied  load  of  66  kN 
(30  MPa)  as  schematically  illustrated  in  Fig.  3.37.  There  was  no  load  drop  when  this 
happened,  and  the  external  wing  cracks  remained  closed. 

A  schematic  illustration  of  the  fracturing  observed  between  60'- 105°  flaws  (true  for 
both  specimens)  at  an  applied  load  of  74  kN  (33  MPa)  is  shown  in  Fig.  3.38.  The  upper 
internal  wing  crack  propagated  towards  the  lower  internal  flaw  tip  but  it  was  not  apparent 
whether  this  wing  crack  had  actually  connected  with  the  flaw  tip.  No  load  drop  occured 
and  the  external  wing  cracks  stayed  closed.  One  specimen  was  unloaded  after  the  fracture 
pattern  in  Fig.  3.38  was  observed.  Cutting  this  specimen  revealed  that  the  wing  crack  had 
not  connected  with  the  lower  internal  flaw  tip,  i.e.,  coalescence  had  not  taken  place.  The 
other  specimen  was  loaded  until  it  failed  explosively  under  an  applied  load  of  84  kN  (37 
MPa). 


The  fracture  coalescence  between  60’- 135'  flaws  is  shown  in  Figs.  3.39(a)  and  (b). 
For  one  specimen  schematically  drawn  in  Fig.  3.39(a),  the  lower  internal  wing  crack 
propagated  towards  the  upper  external  flaw  tip  at  an  applied  load  of  72.5  kN  (32  MPa). 
For  the  other  specimen  shown  in  Fig.  3.39(b),  the  upper  internal  wing  crack  propagated 
towards  the  lower  external  flaw  tip  at  an  applied  load  of  75  kN  (34  MPa).  The  fracturing 
just  described  was  accompanied  by  a  load  drop  in  both  specimens.  Thus,  although  these 
specimens  were  not  cut,  it  is  likely  that  the  flaws  were  connected  by  wing  cracks  because 
of  the  observed  correlation  between  coalescence  and  load  drops  in  the  other  flaw 
geometries. 

Finally,  for  the  60’- 150’  flaws,  no  coalescence  was  observed  for  both  specimens 
which  were  loaded  up  to  about  75  kN  (33  MPa).  At  this  load,  the  upper  internal  wing 
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crack  had  propagated  towards  the  lower  internal  flaw  tip  as  schematically  drawn  in  Fig. 
3.40.  It  was  not  apparent  whether  this  wing  crack  had  actually  connected  with  the  flaw  tip. 
Loading  was  discontinued  at  75  kN  because  it  was  thought  that  any  additional  increase  in 
load  would  have  resulted  in  complete  disintegration  of  the  specimens  such  as  what  occured 
for  one  of  the  specimens  with  60*- 105*  flaws.  This  latter  case  showed  that  no  additional 
information  about  fracture  patterns  could  be  obtained  from  specimens  which  failed 
explosively. 
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Figure  3.33(a)  Internal  wing  crack  growth  resulting  in  coalescence  between  30*-12Q* 
flaws,  specimen  #1  [applied  load  =  62.4kN]. 
int.  wc  =  internal  wing  crack,  ext.  wc  =  external  wing  crack 


Figure  3.33(b)  Internal  wing  crack  growth  resulting  in  coalescence  between  30*- 120* 
flaws,  specimen  #2  [applied  loads  =  61kNl. 
int  wc  =  internal  wing  crack,  ext.  wc  =  external  wing  crack. 
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Figure  3.34  Internal  wing  crack  growth  resulting  in  coalescence  between  30*- 135*  flaws 
specimen  #1,  accompanied  by  spalling  and  cracking  marked  by  "A"  [applied 
load  =  67.5kN], 

int.  wc  =  internal  wing  crack,  ext.  wc  =  external  wing  crack,  sp  =  spalling. 
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"Apparent"  coalescence  due  to  internal  wing  crack  growth  between  30*- 135 
flaws,  specimen  #2  [applied  load  =  53.5kN]. 
int  wc  -  internal  wing  crack. 


Figure  3.36  Secondary  crack  growth  resulting  in  coalescence  between  45*- 135*  flaws, 
specimen  #1  [applied  load  =  72kN].  Note:  secondary  cracks  grow  between 
internal  and  external  flaw  tips. 

int  wc  =  internal  wing  crack,  ext.  wc  =  external  wing  crack, 
sc  =  secondary  crack. 
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Figure  3.37  Wing  crack  growth  from  45*- 135*  flaws,  specimen  #2.  Note:  lower 

internal  wing  crack  propagated  towards  upper  flaw  [applied  load  =  66  kN], 
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3.2.2  Wing  Crack  Initiation  and  Coalescence  Loads 

The  applied  loads  at  which  wing  cracks  became  visible  is  plotted  versus  the 
ligament  inclination  angle  in  Fig.  3.41.  Triangular,  square  and  circular  symbols  were  used 
for  30*,  45*  and  60*  inclined  flaws  respectively.  For  the  same  ligament  inclination  angle, 
the  wing  crack  initiation  loads  were  always  lowest  for  the  30*  inclined  flaws  and  highest 
for  the  60*  inclined  flaws.  An  explanation  for  this  observation,  given  in  more  detail  in 
Section  4.2,  can  be  found  from  the  variation  of  stress  fields  around  flaws  with  different 
inclination  angles.  Furthermore,  for  the  same  flaw  inclination,  there  appears  to  be  some 
dependence  of  the  wing  crack  initiation  loads  on  the  ligament  inclination  angle.  For  all 
three  inclination  angles,  the  wing  crack  initiation  loads  were  lowest  for  ligament  inclination 
angles  between  50°  and  100°. 

Figs.  3.42,  3.43  and  3.44  are  plots  of  maximum  test  loads3-5  versus  ligament 
inclination  angles  for  30°,  45*  and  60*  inclined  flaws,  respectively.  Dotted  lines  divide  the 
graphs  into  "non-overlapping"  (ligament  inclination  angle  £  90’)  and  "overlapping" 
(ligament  inclination  angle  >  90’)  regions.  Hollow  symbols  are  used  to  plot  maximum 
loads  at  which  coalescence  actually  took  place.  As  noted  in  the  previous  section, 
coalescence  always  took  place  through  secondary  crack  growth  for  non-overlapping  flaws, 
and  through  either  wing  or  secondary  crack  growth  for  the  overlapping  flaws.  Solid 


3  5  The  maximum  test  load  is  the  largest  load  applied  during  each  test  In  cases  where  loading  was  stopped 
after  coalescence  (see  Figs.  3.9, 3.13  and  3.19),  the  maximum  est  load  corresponds  to  the  coalescence  load. 
Otherwise,  the  maximum  test  load  is  the  load  at  which  a  specimen  failed  (i.e.,  severely  cracked)  or  at  which 
testing  was  stopped  before  coalescence  or  failure  occured. 
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symbols  are  used  to  plot  maximum  loads  that  correspond  either  to  loads  at  failure  without 
coalescence,  or  loads  at  which  testing  was  discontinued  before  failure. 

In  general,  for  the  same  flaw  inclination  angles,  the  coalescence  loads  were  lower 
for  non-overlapping  when  compared  to  overlapping  flaw  geometries.  The  rest  of  this 
report  will  focus  on  modelling  coalescence  mechanisms  between  non-overlapping  flaws 
since  these  cases  are  usually  more  critical  and  exhibit  secondary  crack  growth  that  cannot 
be  explained  by  existing  fracture  mechanics  theories  (see  Section  2.3).  Fig.  3.45  is  a  plot 
of  the  observed  coalescence  loads  versus  ligament  inclination  angle  for  only  non¬ 
overlapping  flaws  with  30°,  45*  and  60*  inclinations.  One  can  see  that  the  measured 
coalescence  loads  are  consistent  between  specimens  with  identical  flaw  geometries  (two 
were  tested  for  each  geometry).  This  established  the  reproducibility  of  the  observed 
coalescence  phenomenon.  Furthermore,  for  the  same  ligament  inclination  angle, 
coalescence  loads  increased  witflCincreasing  flaw  inclination  angle.  Finally,  for  a  given 
flaw  inclination  angle,  coalescence  loads  were  lowest  for  a  ligament  orientation  angle  of 
75*.  It  was  the  goal  of  the  modelling  effort  described  in  Chapter  4  to  reproduce  these 
trends. 
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Figure  3.41  Wing  crack  initiation  loads. 
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Figure  3.42  Maximum  test  load  vs.  ligament  inclination  angle  for  30*  flaws. 
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Figure  3.43  Maximum  test  load  vs.  ligament  inclination  angle  for  45’  flaws. 
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Figure  3.44  Maximum  test  load  vs.  ligament  inclination  angle  for  60‘  flaws. 
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Figure  3.45  Coalescence  load  vs.  ligament  inclination  angle  for  non-overlapping  flaw 
geometries. 
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3.3.3  Fracture  Surface  Morphologies 


Fig.  3.46  contains  schematic  drawings  of  the  specimens  with  30‘-45*  flaws3-6  after 
they  had  been  cut  to  expose  the  wing  and  coalescence  crack  surfaces.  The  three  types  of 
areas,  distinguished  by  the  markings  on  the  figures,  may  be  described  as  follows: 

Type  I  (dotted  region’)  :  This  is  characterized  by  the  complete  absence  of  pulverized 
gypsum  and  a  smooth  appearance.  A  magnified  view  of  this  surface  type  (Fig.  3.47) 
shows  that  the  surface  is  not  glassy  smooth;  however,  the  uniformity  of  the  roughness 
makes  it  appear  smooth  when  the  surface  is  observed  without  any  magnification. 

Type  II  (dashed  region)  :  This  is  characterized  by  traces  of  pulverized  gypsum.  A 
magnified  view  (see  Fig.  3.48)  shows  that  the  surface  is  not  uniformly  rough  like  the  Type 
I  surface.  Instead,  small  patches  of  roughness  are  interspersed  in  a  surface  area  that  is 
actually  smoother  than  the  Type  I  surface.  However,  the  patchy  roughness  seen  under 
magnification  causes  this  type  of  fracture  surface  to  appear  rougher,  when  viewed  without 
magnification,  than  the  Type  I  surface. 

Type  III  (shaded  region)  :  This  is  characterized  by  significant  amounts  of  pulverized 
gypsum.  A  magnified  view  of  this  surface  is  shown  in  Fig.  3.49.  The  roughness  is 
interspersed  as  in  the  Type  II  surface;  however,  the  ridges  are  higher  in  the  Type  HI  surface 
which  causes  it  too  look  rougher  than  the  Type  II  surface  when  viewed  without 
magnification. 

The  Type  I  and  Type  II  surfaces  were  also  observed  under  a  scanning  electron 
microscope  at  a  ~1000x  magnification  (Figs.  3.50  and  3.51  respectively).  Fig.  3.50 


3  6  Thes  are  the  same  specimens  shown  in  Figs.  3.1 1  and  3.23. 
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shows  the  porous  microstructure  of  gypsum  which  consists  of  interlocking  needle-like 
grains.  When  Figs.  3.50  and  3.51  are  compared,  the  Type  II  surface  shows  more 
compacted  grains  with  less  pore-space  than  the  Type  I  surface.  This  is  probably  why  the 
area  between  the  patchy  roughness  of  the  Type  II  surface  seen  in  Fig.  3.48  appears  to  be 
smoother  than  the  Type  I  surface  in  Fig.  3.47. 

A  photograph  at  ~1000x  magnification,  shown  in  Fig.  3.52,  was  taken  very  close 
to  one  of  the  internal  flaw  tips  where  there  was  a  very  small  region  of  pulverized  gypsum. 
The  micrograph  shows  severely  crushed  and  compacted  grains,  especially  when  compared 
to  the  micrographs  in  Figs.  3.50  and  3.51. 

Fig.  3.53  and  3.54  show  the  surface  features  of  cracks  growing  from  30° -75“  and 
45“ -60"  flaws.  It  is  evident  from  these  figures  that  coalescence  crack  surfaces  are 
consistent  in  blocks  with  identical  flaw  geometries,  as  in  the  case  of  30“-45“  flaws  shown 
in  Fig.  3.46.  Furthermore,  in  comparing  coalescence  crack  surfaces  among  the  different 
flaw  arrangements,  one  sees  that  coalescence  crack  surface  features  vary  with  the  flaw 
arrangement. 

All  the  external  wing  cracks,  regardless  of  the  flaw  inclination  or  the  flaw 
geometry,  only  had  Type  I  surfaces3-7.  Hackle  marks,  which  are  usually  associated  with 
dynamic  crack  propagation,  were  seen  on  the  external  wing  cracks  but  at  some  distance 
from  the  flaw  tips.  It  is  not  surprising  that  these  marks  were  not  present  near  the  flaw  tips 
since  the  wing  cracks  grew  quasi-statically  initially  and  only  extended  rapidly  at 


3,7  It  is  likely  that  the  internal  wing  cracks  had  the  same  surface  type.  However,  the  surfaces  of  these 
cracks  were  not  exposed  by  the  cutting  process  shown  in  Fig.  3.29  for  cases  where  coalescence  was  through 
secondary  cracking. 
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coalescence.  In  contrast  to  the  wing  cracks,  the  coalescence  cracks  consisted  of  regions 
with  different  surface  types,  as  can  be  seen  in  Fig.  3.46.  3.53  and  3.54. 

The  Type  I  surface  is  characteristic  of  surfaces  that  run  parallel  or  subparallel  to  the 
applied  uniaxial  compressive  load.  Examples  of  these  are  the  dotted  regions  in  the  30*-75* 
flaws.  Thus,  the  Type  I  surface  is  typically  found  in  coalescence  cracks  in  specimens  with 
large  ligament  inclination  angles  (i.e.,  close  to  90’).  The  Type  II  surfaces  are  inclined  at 
intermediate  angles  with  respect  to  the  horizontal  while  the  Type  III  surfaces  are  sub¬ 
horizontal  or  are  inclined  at  very  low  angles  with  respect  to  the  horizontal.  Both  these  latter 
types  are  present  in  all  coalescence  cracks  but  with  varying  areal  extents  among  the  flaw 
geometries.  They  are  usually  found  in  coalescence  cracks  in  flaw  geometries  with  smaller 
ligament  orientation  angles  (e.g.,  30’ -45*  and  45*-60’  flaws). 

The  complete  absence  of  pulverized  gypsum  in  the  Type  I  crack  surface  could  be 
considered  as  evidence  that  this  portion  of  the  crack  was  created  in  tension  since  the 
surfaces  must  have  immediately  pulled  apart  after  the  crack  was  formed.  On  the  other 
hand,  the  presence  of  pulverized  gypsum  in  the  Type  II  and  Type  III  crack  surfaces  does 
not  necessarily  mean  that  these  portions  of  the  coalescence  cracks  were  created  in  shear:  it 
is  also  possible  that  they  were  created  in  tension  but  that  the  surfaces  slid  past  each  other 
causing  the  gypsum  to  be  pulverized. 
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Figure  3.46  Surface  morphologies  of  cracks  in  cut  specimens  with  30’-45*  flaws. 


Figure  3.47  Micrograph  (~15  x)  of  Type  I  surface. 


Figure  3.48  Micrograph  (~15  x)  of  Type  II  surface. 


Figure  3 .49  Micrograph  (~  1 5  x)  of  Type  in  surface. 
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Figure  3.50  Micrograph  (1000  x)  of  Type  I  surface. 
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Figure  3.5 1  Micrograph  (1000  x)  of  Type  II  surface. 
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Figure  3.53  Surface  morphologies  of  cracks  in  cut  specimens  with  30'-75°  flaws. 
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Figure  3.54  Surtace  morphologies  of  er.uk s  in  cut  specimens  with  45°-60°  flaw 
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Chapter  4 

Analysis  and  Coalescence  Model  Development 


4.1  Introductory  Remarks 

The  modelling  process  described  in  this  chapter  consisted  of  three  stages.  In  the 
first  stage,  stress  fields  around  interacting  flaws  were  established  using  analytical  elastic 
solutions  for  plates  with  elliptical  holes  under  far-field  uniform  loading  (Section  4.2  i. 
Correlations  were  made  between  these  stress  fields  and  the  observed  loads  at  which  wing 
cracks  initiated.  However,  it  was  noted  in  the  previous  chapter  that  coalescence  always 
came  after  wing  cracks  had  propagated  Thus,  in  the  second  stage  af  the  modelling 
process,  the  stress/strain  fields  around  flaws  with  wing  cracks  were  determined  (Section 
4.3).  The  finite  element  method  was  used  in  the  stress  analyses  since  analytical  solutions 
do  not  exist  for  open  flaws  with  w  ing  cracks.  By  correlating  the  observed  coalescence 
mechanisms  in  the  experiments  with  the  calculated  stress/strain  fields  around  interacting 
flaws  with  wing  cracks,  the  suitability  of  some  fracture  modelling  approaches  was 
evaluated.  A  smeared-crack  damage  mode!  coupled  with  a  strain-based  failure  criterion 
was  found  to  be  most  appropriate  for  modelling  coalescence.  Therefore,  in  the  final  stage 
of  the  modelling  process,  a  strain-based  damage  model  was  formulated,  implemented  and 
used  to  simulate  the  experiments  described  in  Chapter  3  (Section  4.4).  Very  satisfactory 
agreement  between  the  observed  and  simulated  fracture  behavior  w  as  achieved. 
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4.2  Interacting  Elliptical  Flaws  in  an  Infinite  Medium  under 
Far-field  Compression 

4.2.1  Solution  for  a  single  flaw 


The  stress  field  around  an  elliptical  flaw  under  far-field  uniaxial  stress  has  been 
derived  by  Inglis  (1913)  and  is  given  in  elliptical  coordinates  by  the  following  expressions 
(K.  Chang,  1981): 

+  °Tn  =  pe2’c  cos2(3 


+  ap 


sinh 24(l  -e2'"  cos2pj-e2^0  sin 2(3sin 2r| 


(4.1a) 


o,:  -  oT)  =  o  pl-sinh24(cosh2^()  -cos2p)  +  sinh2^cos2r|^l  -e2’“  cos2p) 

-  e2’0  sin2Psin2r)cosh24-e2:i°  sinh2^sinh  2(t,  -  ^0)cos2(ri  -  P) 

-  e2"0  sin2r)cosh2(4-^„)sin2(ri-P)]  (4.1b) 

+  2c tpe2l°  cosh2(^-^0)cos2(r]-P) 


=  —  a:/)[(l  -  e2'~°  cos2p)sin2rjcosh2^-sin  2ri(cosh2^(>  -cos2p) 

+  e2’1  sin2Pcos2r)sinh2c,  +  e2;°  sinh2qcosh2(4-£0)sin2(ri-P) 

„  (4.1c) 

-  e2,c  sin2risinh2(4-s0)cos2(ri-p)] 

-  ape:'°  sin  2(rj  -  P)sinh  2(4  -  £0) 

where  a  =  (cosh  2c  -  cos  2r ]Y\  \  and  rj  are  elliptical  coordinates  as  illustrated  in  Fig.  4.1, 
is  the  elliptical  coordinate  corresponding  to  the  boundary  of  the  elliptical  flaw,  p  is  the 
far-field  stress  (see  Fig.  4.1,  positive  for  tension,  negative  for  compression)  and  P  is  the 
angle  between  the  applied  stress  and  the  major  axis  of  the  elliptical  flaw  (i.e.  it  is  equal  to 
90°  minus  the  flaw  inclination  angle;  the  latter  as  defined  in  Chapter  3). 
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The  observed  trends  in  wing  crack  initiation  loads  (Section  3.3.2)  can  be  compared 
to  predictions  using  the  single  flaw  solutions.  At  the  edge  of  the  flaw,  c =  0  so  that 
can  be  obtained  by  substituting  ^=40  into  Eqn.  (4.1a): 

)  =  a/>[cos2p  +  sinh  2^0  -  e2‘°  cos(p  -  T))].  (4.2) 

This  stress  has  a  maximum  value  at  some  rj  as  shown  in  Fig.  4.2  where  it  is  plotted  against 
r\  for  Hass'  inclination  angles  30°.  45°  and  60°  (corresponding  to  (5  =  60°,  45°  and  30°, 
respectively).  Assuming  that  fracturing  initiates  when  the  maximum  tensile  hoop  stress 
reaches  a  critical  value  (the  maximum  tensile  principal  stress  criterion  or  MTPS)  and  using 
Fig.  4.2,  one  would  predict  that  crack  initiation  will  occur  at  a  lower  load  for  the  60° 
inclined  flaw  when  compared  to  the  45°  and  30°  inclined  flaws.  However,  quite  the 
opposite  was  observed  in  the  experiments  described  in  the  previous  chapter.  The  wing 
cracks  for  the  30°  inclined  flaws  were  visible  at  low-er  loads  when  compared  to  wing 
cracks  growing  from  60°  inclined  flaws  (see  Fig.  3.41  for  wing  crack  inititiation  loads 
from  experiments).  A  possible  explanation  for  this  apparent  disagreement  between  the 
theoretical  prediction  and  the  experimental  observation  will  be  given  below. 

Figs.  4.3(a),  (b)  and  (c)  are  tensile  principal  stress  vector  plots  for  the  three  flaw 
inclinations  30  .  45°  and  60°.  All  three  cases  were  computed  for  the  same  compressive 
load,  the  magnitude  of  which  is  shown  at  the  bottom  of  the  Figures.  These  figures  show 
that  high  tensile  stresses  exist  around  the  flaw  tips  but  eventually  decrease  away  from  the 
flaw  tips.  Furthermore,  even  though  the  hoop  stresses  are  higher  for  the  60°  flaw,  the 
tensile  stresses  decrease  more  rapidly  with  distance  from  the  flaw  tip  for  the  60°  flaw  when 
compared  to  the  45'  and  30°  flaws.  The  wing  cracks  growing  from  the  60°  flaws  may 
initiate  at  lower  loads  but  need  only  propagate  a  short  distance  to  reach  an  area  with  lower 
tensile  stresses  On  the  other  hand,  wing  cracks  growing  from  the  45°  and  30°  flaws 
would  have  to  propagate  a  longer  distance  before  reaching  an  area  with  lower  tensile 
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stresses.  Thus,  it  is  no  longer  surprising  that  wing  cracks  were  visible  at  lower  loads  for 
the  4 5°  and  30c  flaws  even  though  Eqn.  (4.2)  predict  lower  wing  crack  initiation  loads  for 
60°  flaw  s.  This  highlights  the  fact  that  the  determination  of  the  "critical  flaw  inclination"  in 
compressive  fields  must  be  based  on  studies  of  crack  growth  rather  than  crack  initiation 
alone,  as  Hoek  and  Bieniawski  (1965)  and  others  have  pointed  out.  This  is  the  primary 
difference  between  fractjre  studies  in  compression  and  tension:  "failure"  in  tension 
generally  coincides  w  ith  crack  initiation  due  to  unstable  crack  propagation  that  is  frequently 
associated  with  predominantly  tensile  stress  fields.  On  the  other  hand,  propagation  in 
compressive  fields  is  usually  stable,  i.e.,  a  load  increase  is  required  for  cracks  to  lengthen. 


4.2.2  Superposition  of  Stress  Fields  from  Two  Elliptical  Flaws 

A  first-order  approximation  of  the  stress  field  around  two  interacting  elliptical  flaws 
can  be  obtained  by  superposing  the  fields  given  by  Eqns.  (4.1a-c).  This  superposition 
procedure  is  applied  to  various  flaw  geometries  and  the  results  for  30°-75°  and  3()°-45° 
arrangements  are  shown  in  Figs.  4.4  and  4.5,  respectively.  These  are  major  principal 
stress  vector  plots  and  only  the  tensile  stresses  are  plotted.  Note  that  tensile  stresses  exist 
w  ithin  the  ligament  area  between  the  flaw  s  and,  if  a  fracture  criterion  based  on  maximum 
tensile  principal  stresses  is  adopted,  crack  initiation  and  growth  in  that  region  may  be 
expected.  However,  tensile  stresses  of  about  the  same  magnitude  exist  in  areas  outside  of 
the  ligament  (encircled  areas  in  Figs.  4.4  and  4.5)  although  no  such  cracking  was  observed 
in  these  regions  for  30°-45°  (see  photograph  in  Fig.  3.12)  and  30'-75°  flaws.  It  is  possible 
that  these  latter  tensile  stresses  were  relieved  by  wing  crack  growth  which  always  preceded 
coalescence  (see  Section  3.2.1).  However,  the  relief  of  tensile  stresses  within  the  ligament 
by  wing  crack  growth  is  just  as  likely.  Thus,  there  is  obviously  a  need  to  determine  stress 
fields  around  interacting  flaws  after  w-ing  crack  growth  before  any  definite  conclusions  can 
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be  drawn  regarding  the  mechanism  underlying  the  experimentally  observed  fracture 
coalescence.  In  recognition  of  this  need,  a  series  of  finite  element  studies  were  performed 
to  determine  the  stress  fields  around  pre-existing  flaws  with  wing  cracks.  The  results  of 
these  analyses  are  described  in  the  following  section. 
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Figure  4.1  (a 


Figure  4. 1 
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)  Inclined  elliptical  flaw  subjected  to  uniaxial  compression 


)  Definition  of  elliptical  coordinates 
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— —  Scale:  Applied  Vertical  Compressive  Stress 


Figure  4.3(a)  Vector  plot  of  the  maximum  principal  stresses  around  a  30'  flaw  under  far- 
field  uniaxial  compression  (only  tensile  stresses  are  plotted). 


"“■*  Scale:  Applied  Vertical  Compressive  Stress 


Figure  4.3(b)  Vector  plot  of  the  maximum  principal  stresses  around  a  45*  flaw  under  far- 
field  uniaxial  compression  (only  tensile  stresses  are  plotted). 
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Scale:  Applied  Vertical  Compressive  Stress 

Figure  4.3(c)  Vector  plot  of  the  maximum  principal  stresses  around  a  60'  flaw  under  far- 
field  uniaxial  compression  (only  tensile  stresses  are  plotted). 
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Scale  Applied  Vertical  Compressive  Stress 


Figure  4.4 


Vector  plot  of  maximum  principal  stresses  around  30’-75"  flaws  under  far- 
field  uniaxial  compression  (obtained  from  superposing  single  flaw  solution; 
only  tensile  stresses  are  plotted;.  High  tensile  stresses  away  from  ligament 
area  exist  within  encircled  region. 
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Scale:  Applied  Vertical  Compressive  Stress 


Figure  4.5 


Vector  plot  of  maximum  principal  stresses  around  30'-45*  flaws  under  far- 
field  uniaxial  compression  (obtained  from  superposing  single  flaw  solution; 
only  tensile  stresses  are  plotted).  High  tensile  stresses  away  from  ligament 
area  exist  within  encircled  region. 
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4.3  Stress  Analysis  of  Interacting  Flaws  with  Wing  Cracks 

4.3.1  Description  of  Finite  Element  Models 

In  order  to  establish  the  stress  conditions  after  wing  crack  growth  and  prior  to  crack 
coalescence,  finite  element  analyses  were  performed  on  a  number  of  the  flaw  arrangements 
that  were  studied  experimentally.  A  linear  elastic  material  was  assumed  for  these  analyses, 
primarily  because  there  were  no  signs  of  gross  material  plasticity  or  damage  prior  to 
coalescence  in  the  experiments.  A  Young's  modulus  of  approximately  6200  MPa  was 
obtained  from  the  experimental  load-displacement  curves  and  used  in  the  analyses.  A 
Poisson's  ratio  of  0.28  was  assumed. 

Fig.  4.6  shows  an  example  of  a  finite  element  model  for  determining  the  stress  field 
around  pre-existing  flaws  (in  this  case.  30e-75‘  flaws)  with  wing  cracks.  It  must  noted  that 
no  attempt  was  made  to  simulate  the  actual  growth  of  the  wing  cracks.  Successful 
simulations  of  wing  crack  growth  have  already  been  reported  by  Ingraffea  and  Heuze 
( 19S0).  Furthermore,  the  experiments  in  Chapter  3  have  show  n  that  the  growth  of  wing 
cracks  did  not  directly  result  in  coalescence  although  it  may  have  indirectly  affected  the 
crack  growth  by  altering  the  stress  field  around  the  pre-existing  flaws.  The  elements  on 
each  side  of  the  thickened  lines  in  the  finite  element  model  shown  in  Fig.  4.6  are  not 
"connected"  by  common  nodes,  thus  forming  a  discontinuity  corresponding  to  the  wing 
cracks  along  these  lines.  The  shape  of  the  wing  cracks  in  the  analyses  were  obtained  from 
photographs  of  wing  cracks  taken  during  the  experiments. 

All  elements  in  the  model  are  plane  strain,  8-noded  quadrilaterals.  Special  singular 
elements  were  not  used  around  the  wing  cracks  tips  because  the  stress  intensity  factors  at 
the  wing  crack  tips  were  only  of  secondary  importance  in  this  study. 
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The  following  boundary  conditions  were  applied  on  the  finite  element  model:  zero 
vertical  displacements  along  the  bottom  edge,  and  a  uniform  distributed  load  on  the  top 
surface,  the  magnitude  of  which  was  approximately  the  measured  coalescence  load  (in 
MPa)  for  the  flaw  geometry  being  analyzed  (see  Section  3.3.2).  It  was  assumed  that  the 
flaws  were  sufficiently  far  from  the  top  and  bottom  edges  of  the  block  such  that  the  exact 
distribution  of  loads  on  these  edges  did  not  significantly  affect  the  stresses  around  the 
flaws.  Because  the  model  is  two-dimensional,  three-dimensional  effects  in  the  experiments 
were  neglected. 

4.3.2  Linear  Elastic  Analysis  Results  and  Discussion 

The  analyses  showed  that  tensile  stresses  in  the  ligament  areas  were  not  relieved  by 
wing  crack  growth,  a  possibility  that  was  raised  in  Section  4.2.2.  This  can  be  seen  in 
Figs.  4.7  and  4.N  w  hich  are  major  principal  stress  vector  plots  generated  from  the  analysis 
results  for  3(>'-75'  and  30‘-45*  flaws  (compare  with  Fig.  4.4  and  4.5). 

The  stress  fields  after  wing  crack  growth  were  combined  with  experimental 
observations  to  look  into  possible  models  that  would  reproduce  the  observed  coalescence 
process.  1  wo  possible  approaches  to  modelling  the  coalescence  fracture  are:  (1 )  a  discrete 
crack  approach  coupled  with  the  maximum  tensile  principal  stress  fracture  criterion 
(MTPS).  and  (2)  a  smeared  crack  approach  coupled  with  MTPS.  The  appropriateness  of 
each  approach  will  now  be  discussed 

In  a  discrete  crack  approach  coupled  with  MTPS.  cracks  are  assumed  to  propagate 
perpendicular  to  the  direction  of  maximum  tensile  principal  stress.  The  coalescence  cracks 
shown  in  the  photographs  of  30°-75°  flaws  (see  Figs.  3.24  and  3.25)  do  lie  perpendicular 
to  the  stress  vectors  plotted  in  Fig.  4.7.  Thus,  in  the  case  of  the  30"-75’  flaws,  the  discrete 
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crack  approach  coupled  with  MTPS  may  succesfully  simulate  coalescence.  However,  the 
path  of  the  coalescence  cracks  shown  in  the  photographs  of  30°-45°  flaws  (see  Figs.  3.23 
and  3.27 )  do  not  run  perpendicular  to  the  tensile  stress  vectors  in  Fig.  4.8.  This  is  contrary 
to  what  is  assumed  when  applying  the  discrete  crack  approach. 

In  a  smeared  crack  approach  coupled  with  MTPS,  cracks  initiate  and  propagate 
along  regions  where  the  highest  tensile  principal  stresses  exist.  However,  in  contrast  to  the 
discrete  crack  approach,  these  cracks  do  not  necessarily  grow  perpendicular  to  the  direction 
of  the  tensile  stresses.  If  the  magnitude  of  the  maximum  tensile  principal  stress  is  used  as  a 
criterion  for  smeared  cracking,  this  approach  would  predict  cracks  initiating  within  the 
shaded  regions  on  the  stress  contour  plots  for  the  30“ -75'  and  30“-45°  flaws  (see  Figs.  4.9 
and  4.10).  The  approximate  path  of  the  observed  coalescence  crack  betw  een  the  30'-75° 
flaws,  superposed  on  the  stress  contour  plot  in  Fig.  4.9,  does  intersect  this  shaded  region. 
On  the  other  hand,  this  is  not  true  for  coalescence  crack  between  the  30° -45'  flaws  which  is 
schematically  show  n  in  Fig.  4.10.  Further  evidence  against  the  use  of  the  smeared  crack 
approach  coupled  with  MTPS  is  shown  in  Fig.  4.1 1(a).  This  is  an  expanded  view  of  the 
stress  field  around  the  30c-45°  flaws  in  Fig.  4.10.  Regions  of  high  tensile  stresses  exist 
outside  of  the  ligament  between  the  flaws  'refer  to  dotted  areas  in  Fig.  4.1 1(a)).  Thus,  a 
smeared  crack  approach  coupled  with  MTPS  would  predict  cracking  in  these  areas  whereas 
no  such  cracking  was  observed  in  the  experiments  (see  photograph  in  Fig.  3.12).  Unless 
the  actual  stress  field  is  different  from  the  computed  stress  Field  in  the  case  of  the  3()'-45° 
flaws,  it  appears  that  the  smeared  crack  approach  using  MTPS  will  not  successfully 
simulate  the  coalescence  crack  for  this  geometry. 

A  third  possibility  for  modelling  fracture  coalescence  is  a  smeared  crack  approach 
coupled  with  a  maximum  tensile  principal  strain  fracture  criterion  (MTPSN).  The 
suitability  of  this  approach  can  also  be  tested  against  the  observed  fracture  patterns  and  the 
results  of  the  stress  analyses  described  above.  Fig.  4.1 1(b)  is  a  contour  plot  of  the  major 


principal  strains  from  the  same  analysis  that  was  used  to  generate  the  principal  stress 
contour  plot  in  Fig.  4.11(a)1.  Aside  from  areas  around  the  wing  crack  tip  where  a 
singularity  was  introduced  by  the  sharpness  of  the  crack  tip,  regions  of  highest  tensile 
strains  exist  only  at  the  tips  of  the  pre-existing  flaws  and  within  the  ligament  zone.  Thus, 
a  smeared  crack  approach  based  on  MTPSN  will  predict  cracking  to  initiate  at  all  the  flaw 
tips  and,  from  the  internal  flaw  tips,  to  propagate  through  the  ligament  zone  forming  a 
coalescence  crack  (see  dotted  regions  in  Fig.  4.11(b)).  In  the  experiments,  it  was 
impossible  to  verify  that  the  coalescence  cracks  initiated  at  the  flaw  tips  due  to  the  rapidity 
with  which  these  appeared  in  most  flaw  geometries,  unless  the  observed  spalling  is  taken 
as  an  indicator  of  smeared  cracking  at  the  flaw-  tips.  Flowever,  the  observed  coalescence 
crack  does  approximately  fall  within  the  highly  strained  (dotted)  region  marked  in  Fig. 
4.11(b). 

It  is  also  interesting  to  compare  the  spreading  of  cracked  regions  predicted  by 
tensile  stress-based  and  tensile  strain-based  smeared  crack  models.  Without  running  a  non¬ 
linear  analysis,  one  could  obtain  an  idea  of  this  spreading  from  the  linear  elastic  stress  and 
strain  contour  plots  as  shown  in  Figs.  4.12  through  4.17.  In  these  plots,  regions  between 
stress/strain  contour  lines  that  would  "connect"  the  internal  flaw  tips  have  been  shaded. 
These  shaded  regions  would  then  approximate  the  extent  of  smeared  cracking  prior  to 
coalescence.  In  all  the  flaw  geometries,  the  estimated  smeared  crack  regions  in  the 
principal  strain  contour  plots  (Figs.  4.12(a)  through  4.17(a))  are  more  localized  than  those 
in  the  principal  stress  contour  plots  (Figs.  4.12(b)  through  4.17(b)).  Thus,  the  strain¬ 


'll  must  be  emphasized  that  the  maximum  principal  strain  is  not  simply  proportional  to  the  maximum 
principal  stress  in  the  same  direction,  c.g.,  t\  =  ~  -  jc  (02  +  03).  This  implies  that  tensile  principal 

strums  arc  higher  for  regions  where  there  is  a  combination  of  large  tensile  and  compressive  stresses. 
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bused  smeared  crack  approach  would  simulate  a  more  localized  damaged  region  that  is 
closer  to  w  hat  was  observed  in  the  experiments.  A  similar  observation  regarding  more 
localized  predictions  of  damage  for  a  strain-based  criterion  was  made  by  Tvergaard  (1982) 
in  his  study  of  void  coalescence  in  ductile  materials. 
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Figure  4.6  Finite  element  model  for  determining  stress  field  around  3CT-75’  flaws  w  ith 
wing  cracks. 
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YH 


Scale:  Applied  Vertical  Compressive  Stress 


Figure  4.7  Vector  plot  of  maximum  principal  stresses  around  30"-75*  flaws  with  wing 
cracks  under  far-field  uniaxial  compression  (obtained  from  Finite  element 
analysis  results;  only  tensile  stresses  are  plotted). 


136 


Scale:  Applied  Vertical  Compressive  Stress 


Figure  4.8  Vector  plot  of  maximum  principal  stresses  around  30'-45'  flaws  with  wing 
cracks  under  far-field  uniaxial  compression  (obtained  from  finite  element 
analysis  results;  only  tensile  stresses  are  plotted). 
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4.3.3  A  Comparison  Between  Stress  and  Strain-Based  Failure  Criteria 

One  interesting  implication  of  using  strains  rather  than  stresses  is  the  "automatic" 
incorporation  of  the  biaxial  strength  effect  which,  in  the  case  of  combined 
tensile/compressive  stress  states,  refers  to  a  decrease  in  the  tensile  stress  at  failure  due  to  a 
compressive  stress  that  is  applied  normal  to  the  direction  of  tension.  According  to  the 
MTPS  fracture  criterion,  fracturing  occurs  at  a  point  where  the  maximum  principal  stress 
reaches  a  critical  value,  say  ot.  On  the  other  hand,  the  MTPSN  fracture  criterion  assumes 
that  fracturing  occurs  at  a  point  where  the  maximum  principal  strain  reaches  a  critical  value, 
say  et.  Both  of  these  values  can  be  obtained  from  a  uniaxial  tensile  test  such  that,  if  the 
material  is  elastic  prior  to  fracture,  the  critical  strain  and  stress  are  simply  related  by  = 
EC[  where  E  is  ’he  modulus  of  elasticity. 

Consider  now  the  biaxial  stress  condition: 

o3  =  ~kc, 

where  k  is  a  positive  number.  According  to  MTPS,  fracturing  will  occur  when  C]  =  ot. 
regardless  of  the  value  of  03.  However,  a  different  result  will  be  obtained  with  MTPSN 
since  both  Oj  and  03  will  be  contributing  to  the  maximum  principal  strain  as  follows: 

el=E  (°|-v03)=E(1  +  Vk)o.  (4.3) 

but  since  ot  =  Ect.  MTPSN  will  predict  fracturing  at 


Since  k  is  a  positive  number,  then  the  fracture  stress  according  to  MTPSN  is  lower  than  the 
fracture  stress  according  to  MTPS.  Furthermore,  Eqn.  (4.4)  implies  that  the  fracture  stress 
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is  reduced  when  k  is  not  equal  to  zero.  Thus,  the  biaxial  strength  effect  follows  from  using 
the  maximum  principal  strain  as  a  failure  criterion. 

If  both  stresses  are  tensile,  i.e.,  k  is  a  negative  number,  then  Oj  would  be  greater 
than  in  Eqn.  (4.4).  Thus,  the  use  of  MTPSN’  for  biaxial  stress  conditions  where  both 
stresses  are  tensile  would  result  in  a  fracture  stress  that  is  higher  than  the  uniaxial  fracture 
stress.  However,  no  experimental  data  on  the  biaxial  tensile  strength  of  rocks  exist  to 
support  or  refute  this  result  Fortunately,  in  this  study,  there  are  virtually  no  such  biaxial 
tensile  stress  conditions  around  the  flaws,  except  possibly  around  the  wing  crack  tips 
where  the  stresses  are  singular.  Within  the  ligament  areas  of  the  non-overlapping  flaws, 
the  stress  conditions  are  such  that  always  one  of  the  principal  stresses  is  compressive. 


Another  interesting  consequence  of  using  MTPSN  is  the  possibility  of  tensile 
fracture  in  completely  compressive  Fields.  Consider  an  axisymmetric  stress  state  given  by: 


(4.5) 


where  p  and  q  are  positive  values.  Assuming  that  the  material  is  elastic,  the  maximum 
principal  strain  is: 

E;  =  e2  =  £  (<*:  -  va:  "  VC3)  =  ^  [q(v  -  1)  +  Vp],  (4.6) 


Thus,  the  principal  strain  £]  is  tensile  if: 

p  1-v 


Table  4.1  shows  values  of  minimum  p/q  ratios  for  which  tensile  fracturing  could  occur  in  a 
completely  compressive  stress  field  according  to  MTPSN. 
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Substituting  £]  =  £t  in  Eqn.  (4.6)  and  rearranging,  an  expression  for  the  failure 
envelope  in  (p  -  q)  vs  q  space  is  obtained: 


p  -  q  =  +  ^(1  -  2v) 

v  v 


(4.8) 


If  q  =  0  in  Eqn.  (4.2).  then  p  =  uniaxial  compressive  strength  =  oc: 


P  = 


(4.9) 


Using  Eqn.  (4.9),  a  normalized  failure  envelope  can  be  obtained  from  Eqn.  (4.8): 


p-q  =  ]  +  Ji_ H  ~2v) 

Oc  O,  V 


(4.10) 


The  above  failure  envelope,  with  v  =0.15  and  v=0.28,  is  compared  to  the  Griffith  and 
Modified  Griffith  criteria,  as  well  as  experimental  data  summarized  by  McClintock  and 
Walsh  (1962)  in  Fig.  4.18.  It  is  rather  striking  that  the  failure  envelope  constructed  using 
the  MTPSN  fracture  hypothesis  Fits  the  data  as  well  as  the  Modified  Griffith  criterion. 
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Figure  4. 1 8  Comparison  of  Maximum  Tensile  Principal  Strain,  Modified  Griffith  and 
Griffith  failure  criteria  with  data  from  literature  (from  McClintock  and 
Walsh,  1962).  Note:  a0,  K  and  m  are  material  parameters  required  by  the 
Modified  Griffith  criterion,  oc  is  the  uniaxial  compressive  strength. 
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4.4  A  Tensile  Strain-Based  Damage  Model  for  Coalescence 

It  has  been  shown  in  previous  sections  that  the  smeared  crack  approach  coupled 
with  the  maximum  tensile  principal  strain  criterion  (MTPSN)  seemed  to  provide  a  good 
framework  for  modelling  coalescence.  The  development,  implementation  and  application 
of  a  strain-based  smeared  crack/damage  model  for  coalescence  is  described  below . 


4.4.1  Model  Formulation  and  Implementation 


The  underlying  assumption  of  damage  mechanics  is  that  material  damage  causes  a 
degradation  of  the  material  stiffness.  For  uniaxial  conditions,  it  is  assumed  that  the 
stiffness  degrades  according  to  the  following  rule  [Kachanov,  1958;  Lemaitre.  1986): 

o  =  (l-D)Ee  (4.11) 


where  o  =  uniaxial  tensile  stress,  c  =  uniaxial  tensile  strain,  E  =  elastic  modulus  of  the 
undamaged  material,  and  D  =  a  scalar  internal  state  variable  which  measures  damage.  The 
values  of  D  can  range  from  0  through  1,  corresponding  to  "no  damage"  through 
"completely  damaged"  states.  For  a  strain-based  damage  model,  the  variable  D  is  assumed 
to  be  a  direct  function  of  strain  and.  in  this  study,  the  following  functional  form  for  D  w  as 
used: 

fO  fore  <  E0 

D  =  \  (4  12) 

[H(0  fore  >  eD  K  ’ 


where  H(C)  is  defined  as  follows: 


H(c)  =  1  - 


1 


r 

f  c  M 

exp 

> 

>  \- 

i 

(4.13) 


152 


This  function  has  the  following  mathematical  characteristics:  (1)  when  £  <  e0,  D  =  0,  and 
(2>  as  r  — >  oc  p  1.  The  pammeter  e0  is  a  damage  threshold  value  while  the  parameter 
A,  determines  the  rate  at  which  D  approaches  1  as  strain  increases.  The  parameter  A  is 
limited  to  values  greater  than  0  and  less  than  or  equal  to  l2.  Plots  of  D  vs  £  for  A  =  0.5, 
0.75  and  1.0  are  shown  in  Fig.  4.19  for  £0  =  0.003.  Note  that  larger  A  values  cause  a 
more  rapid  increase  in  damage  with  increasing  strain. 


Substituting  Eqn.  (4.12)  into  Eqn.  (4.13),  an  expression  for  the  variation  of  stress 
versus  strain  is  obtained: 


Ec 


Ec 


(  £  \\ 

exp 

A 

Uo  ) 

for  £  <  £0 
for  £  >  £0 


(4.14) 


This  stress-strain  relation  is  plotted  in  Fig.  4.20  for  £0  =  0.003,  and  A  =  0.5,  0.75,  1.0. 
As  expected,  higher  A  values  give  rise  to  more  pronounced  strain-softening  behavior. 

The  damage  model  defined  above  had  to  be  revised  in  order  to  account  for  the  fact 
that  damage  is  irreversible,  i.e.,  damage  does  not  decrease  with  material  "unloading".  The 
simplest  possible  model  that  accounts  for  damage  irreversibility  is  to  assume  that  the 
material  unloads  along  a  linear  elastic  line  defined  by  the  damage  state  prior  to  unloading. 
To  illustrate  this  loading-unloading  concept,  consider  the  stress-strain  curve  in  Fig.  4.21  on 


2  Justification  for  this  limitation  in  A-valucs  will  be  given  later. 
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which  the  strain  state  at  time  t  ?,  £(t),  denoted  by  £*  is  marked.  Since  £*  >  e0.  the 
corresponding  damage,  D(e*)  can  be  found  from  Eqns.  (4.12)  and  (4.13).  If  the  strain 
deceases  after  time  t,  then  the  materiai  will  foliow  the  linear  curve  given  by  a  =  [3- 
D(£*  )]Ee  (see  Fig.  4.21).  If  the  strain  increases  again  after  unloading  (reloading),  then  the 
material  will  still  follow  the  curve  defined  by  c  =  [1-D(£*)]E£  until  £  becomes  greater  than 
£*.  This  implies  that  damage  is  irreversible  during  unloading  and  that  while  £  <  £*,  no 
additional  damage  occurs  during  reloading.  Therefore.  £*,  the  strain  prior  to  unloading 
which  can  also  be  defined  as  the  maximum  strain  in  the  loading  history  of  the  material,  is 
like  a  "moving”  damage  threshold.  Once  £  >  £*,  damage  starts  to  increase  again  and  the 
material  follows  the  non-linear  curve  as  shown  in  Fig.  4.21. 

The  damage  model  which  accounts  for  these  loading/unloading  conditions  is 
defined  below  [following  Sirno  and  Ju,  1987J: 

1 )  Define  £*(t)  =  a  damage  threshold  that  varies  with  time  t.  Its  initial  value  is  £0. 

2)  Define  a  damage  surface,  analogous  to  the  yield  surface  in  plasticity,  such  that 
strain  states  must  always  be  within  this  surface: 

£(t)-£*(t)<0  (4.15) 

This  implies  that  once  £(t)  exceeds  £0,  £(t)  ->  £*(t),  i.e.,  the  damage  threshold  is  set 
to  the  maximum  strain  in  the  loading  history'  of  the  material. 


3Thc  variable  t  refers  to  "pseudo-time"  in  the  loading  history  [Bathe,  1982],  It  has  nothing  to  do  with 


viscous  effects. 
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3) 


Damage  evolution  is  determined  by  the  following  rule: 


£ 


9H(r) 

df 


D  = 


0 


if  £(t)  =  £  *  (t)  (4.16) 


if  £{ t)  <  E*(t) 


where  1 1(e)  is  defined  by  Eqn.  (4.13).  Eqn.  (4.16)  is  equivalent  to  specifying  that  damage 
increases  (i.e.,  D  >  0)  only  when  the  strain  exceeds  the  current  damage  threshold  £*(t). 
Eqn.  (4.12)  still  holds  but  is  stated  in  differentiated  form  in  Eqn.  (4.16). 

The  uniaxial  constitutive  law  defined  by  Eqn.  (4.11)  can  be  generalized  to 
multiaxial  stress/strain  conditions  by  replacing  the  scalars  by  their  tensorial  representations: 

a,,  =  (1  —  D)Cljk,£k:  (4.17) 

Here.  Cjjiq  are  the  components  of  the  elastic  constitutive  tensor.  Note  that  the  damage 
variable.  D.  is  still  a  scalar.  Thus,  the  model  used  in  this  study  assumes  isotropic  damage, 
as  opposed  to  anisotropic  models  in  the  literature  which  assume  that  damage  can  be 
represented  by  higher  order  tensors  (e.g..  Krajcinovic,  1986;  Valanis,  1990].  The  damage 
evolution  equations  also  need  to  be  generalized  for  multiaxial  conditions.  This  can  be  done 
by  replacing  £  in  Eqn.  (4.16)  with  £cq.  a  scalar  function  of  the  strain  components 
sometimes  referred  to  as  an  equivalent  strain  [Simo  and  Ju,  1987].  Thus,  for  multiaxial 
stress/strain  states,  the  damage  evolution  equations  are: 
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D  - 


an^) 

de_ 


E«,  =E*(t) 


<  £*(t) 


(4.18.) 


where  Hu  )  =  1  - 


exp 

r  Y| 

A  -^-1 

Uo  )  J 

In  this  study,  the  equivalent  strain  was  defined  as  follows: 

ft,  fore,  >  0 


0 


for  e,  <  0 


(4.19) 


where  fi  =  major  principal  strain.  Eqn.  (4.19)  implies  that  only  tensile  strains  cause 
material  damage  in  accordance  with  the  MTPSN  fracture  hypothesis. 

The  preceding  damage  model  was  implemented  as  a  material  subroutine  of  the  finite 

element  program  ABAQUS  [ABAQUS  User's  Manual,  1990],  The  latter  program  uses  a 

Ncwton-Raphson  scheme  to  solve  non-linear  problems.  Therefore,  the  material  subroutine 

must  return  an  incremental  material  stiffness  matrix  to  the  main  program.  By  differentiating 

Eqn.  (4. 17).  an  expression  for  the  material  Jacobian,  — ,  may  be  obtained: 

d£kl 

=  f  1 —  D)Cljk;  —  DC1Jk!  £kJ  (4.20) 


By  using  Eqn.  (4.18)  and  letting 

°fj  =  Cijk!  Ekl 


which  is  the  stress  corresponding  to  an  undamaged  material,  Eqn.  (4.20)  becomes: 
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(4.21) 


o.  =  (1  -  D)Cia.  c.  - 


dH(E«.) 


d£. 


O  ,,E. 


But  since 


Be*,  . 


Be. 


(4.22) 


Eqn.  (4.21)  can  be  written  as 


a  =  (l  -  D)C  .  e. 


3^  ,J  fcu  1 


so  that  the  material  Jacobian  is: 


Bo 


dc 


-  =  (l-D)C,i;- 


dH(e«=) 

Be_  '•  Bck, 


(4.23) 


(4.24) 


Note  that  this  Jacobian  is  not  symmetric  due  to  the  presence  of  tht  second  term. 
Consequently,  there  would  be  a  significant  increase  in  computational  effort  to  solve  the 
resulting  unsymmetric  system  of  equations.  On  the  other  hand,  using  an  approximate  but 
symmetric  Jacobian: 

da,; 

T-1’  O-DJCiju  (4,25) 

dckl 


which  neglects  the  unsymmetric  pan  of  Eqn.  (4.24)  would  lead  to  a  solution  that  would 
require  more  iterations  to  achieve  convergence.  However,  the  increase  in  number  of 
iterations  could  be  offset  by  the  reduced  amount  of  time  required  to  solve  a  symmetric 
system  of  equations.  Thus,  the  approximate  Jacobian  in  Eqn.  (4.25)  was  used  in  the 
implementation  of  the  damage  model.  Furthermore,  Eqn.  (4.25)  essentially  gives  a  secant 
material  stiffeness,  as  opposed  to  a  tangential  material  stiffness  given  in  Eqn.  (4.24). 
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During  initial  attempts  to  use  the  above  material  model  coupled  with  the 
approximate  Jacobian,  difficulties  in  convergence  were  encountered  whenever  a  material 
point  started  to  strain-soften,  i.e.,  unreasonably  small  load  steps  were  required  to  achieve 
convergence.  It  is  possible  that  the  approximate  Jacobian  in  Eqn.  (4.25)  is  not  suitable  for 
solving  strain-softening  problems.  However,  it  is  also  likely  that  the  numerical  schemes 
programmed  into  ABAQUS  do  not  include  a  method  that  is  appropriate  for  material 
softening  [e.g.  indirect  displacement  control  proposed  by  de  Borst,  1986]. 

In  order  to  avoid  the  numerical  problems  associated  with  strain-softening  materials, 
the  following  simplification  was  made.  Rather  than  letting  D  — »  1,  a  saturation  value  of  D 
was  assumed  to  exist  and  correspond  to  the  damage  at  the  peak  of  the  stress-strain  curve. 
By  maximizing  the  second  part  of  Eqn.  (4.14),  one  finds  that  the  peak  stress  occurs  at  £  = 
CsiH  =  f\VA  which  is  given  by  the  following  expression: 


Ef , 


1 


A  exp(l-A) 


(4.26) 


where  the  subscript  "sat"  has  been  used  because  this  also  corresponds  to  the  stress  at  which 
the  saturation  damage  is  reached.  The  reasons  behind  the  previously  stated  constraints  set 
on  possible  values  of  A  are  now  evident.  The  parameter  A  must  be  less  than  1  so  that  £Sat 
>  £0.  Otherwise,  if  A  were  greater  than  1,  a  physically  nonsensical  situation  will  arise 
wherein  damage  saturation  will  occur  at  a  strain  that  is  lower  that  the  strain  at  which 
damage  initiates.  Also  one  can  infer  from  Eqn.  (4.26)  that  A  cannot  be  zero  since  this  will 
produce  a  singular  and  that  A  must  be  a  positive  number  so  that  positive  strains  will 
correspond  to  positive  stresses.  By  substituting  £  =  Esat  =  Eo/A  in  Eqns.  (4.12)  and  (4.13), 
the  saturation  damage  is  obtained: 
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1 


(4.27) 


I 

5 


Dsa.=l- 


exp(l  -  A) 


Therefore,  with  this  assumption,  no  further  damage  occurs  when  £eq  >  (i.e.,  cCq  is 

greater  than  the  strain  at  damage  saturation)  and  the  material  behavior  remains  linearly 
elastic  with  a  reduced  modulus: 

a,,  =  fl  -  DsaJC^i  ey.  (4.28) 

The  uniaxial  stress-strain  curve  for  this  revised,  non-softening  model  is  shown  in  Fig. 
4.22.  Other  authors  have  also  proposed  a  saturation  value  for  damage  for  microcracking 
materials,  with  corresponding  to  a  saturation  microcrack  density  [Charalambides  et  al„ 
1987;  Brockenbrough  et  al„  1987,  Ortiz,  1 987 J. 


Computational  and  uniqueness  problems  associated  with  strain-softening  materials 
were  avoided  in  this  study  through  the  use  of  a  non-softening  damage  model.  Simulations 
of  the  coalescence  processes  observed  in  the  experiments  are  presented  in  the  following 
section. 
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4.4.2  Coalescence  Simulations  Using  a  Strain-based  Damage  Model 

The  objective  of  the  following  simulations  is  to  verify  whether  or  not  the  strain- 
based  damage  model  presented  in  section  4.4.1  can  successful)'  reproduce  the  observed 
coalescence  in  the  experiments  and  the  corresponding  loads  at  which  coalescence  occured. 
The  approach  taken  in  these  simulation  studies  was  as  follows:  (1)  parameters  e0  and  A 
were  estimated  and  modified  until  the  simulated  coalescence  load  was  close  to  the  observed 
value  for  one  flaw  geometry,  in  particular,  for  30°-75°  flaws,  (2)  these  parameters  were 
then  used  to  predict  coalescence  loads  for  the  other  flaw  geometries  and,  (3)  the  predicted 
coalescence  loads  were  then  compared  to  the  observed  values.  Satisfactory  predictions  of 
secondary  crack  growth/coalescence  and  coalescence  loads  as  functions  of  flaw  geometry 
were  achieved. 


The  "Edge"  Finite  Element  Model 

Because  of  limited  computer  resources  and  the  expected  iterative  nature  of  finite 
element  analyses  involving  material  non-linearity,  a  simplified  block  geometry  was  used  in 
applying  the  damage  model  .  "Edge"  models,  where  the  edge  flaws  had  lengths  equal  to 
half  the  flaw  lengths  in  the  full  model,  were  thought  to  be  adequate  for  investigating 
coalescence  crack  growth  in  the  ligaments  (see  Figs.  4.23  and  4.24  for  a  comparison  of 
"edge"  and  "full"  models  of  30° -45°  and  30°-75°  flaw  geometries,  respectively).  At  the 
same  time,  the  "edge”  models  required  a  smaller  number  of  finite  elements  than  the  full 
models.  A  comparison  between  the  strain  fields  using  the  "edge"  and  full  models  are 
presented  in  Figs.  4.25  and  4.26  for  30°-45°  and  30°-75°  flaws  respectively.  Note  that  the 
strain  fields  for  the  30°-45°  flaws  compare  quite  well  whereas  the  strain  levels  for  the  30°- 
75°  flaws  are  higher  in  the  "edge"  model  than  in  the  full  model.  It  was  generally  observed 
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that  edge  models  overestimated  the  strain  levels  in  the  ligaments  when  the  resulting  "edge" 
blocks  were  too  thin,  as  in  the  case  of  the  30°-75“  flaws.  Intuitively,  this  is  not  surprising 
since  the  edges  were  probably  having  a  dominant  effect  on  the  strain  distribution  in  the 
ligament  areas  of  thinner  edge  models. 


Coalescence  Criterion 

A  simulation  of  an  edge  model  using  a  strain-softening  damage  model  would  most 
likely  result  in  a  lcaH  displacement  curve  that  exhibits  a  "peak"  when  the  load-carrying 
capacity  of  the  ligament  has  been  reached.  This  "peak"  load  could  then  be  taken  as  the 
coalescence  load  for  the  edge  model.  However,  no  such  peaks  are  expected  from  the  load- 
displacement  curves  when  a  non-strain-softening  damage  model,  such  as  the  one  proposed 
in  this  study  (see  Fig.  4.22),  is  used.  Therefore,  the  coalescence  criterion  adopted  for  the 
simulations  was  as  follows:  the  coalescence  load  is  the  load  at  which  a  region  of  damage 
with  D  =  Dsj(  connects  the  flaw  tips.  The  application  of  this  criterion  will  be  made  clearer 
in  the  subsequent  discussion  of  simulation  results. 


Parameter  Selection 

In  selecting  the  damage  threshold  value,  e0>  some  guidance  was  obtained  from  the 
linear  elastic  analysis  results  in  Section  4.3.  As  mentioned  in  4.3.1,  strain  fields  in  Figs. 
4.12(a)  through  4.17(a)  were  induced  by  applying  uniformly  distributed  loads  which  were 
approximately  equal  to  the  observed  coalescence  loads.  Thus,  neglecting  any  effects  of 
damage,  these  figures  show  the  approximate  strain  fields  just  prior  to  coalescence.  It  was 
observed  that  the  computed  strains  within  the  ligament  areas  were  within  the  range  of 
0.0022  to  0.0032  from  the  analysis  results  shown  in  Figs.  4.12(a)  through  4.17(a). 
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Therefore,  for  the  first  simulation,  e0  was  set  to  0.003.  A  value  of  0.25  was  arbitrarily 
chosen  for  A.  The  same  modulus  of  6200  MPa  and  Poisson's  ratio  of  0.28  used  in  the 
linear  elastic  analyses  were  adopted  for  the  non-linear  analyses.  Damage  saturation  occurs 
at  a  strain  of  esal  =  eJA  -  0.012  so  that,  from  Eqns.  4.27  and  4.28,  Dsat  =  0.53  and  osat  = 
35  MPa.  The  resulting  stress-strain  curve  is  shown  in  Fig.  4.27. 

Figs.  4.28( a)-(d)  show  the  contours  of  the  damage  variable  D  for  increasing  applied 
loads  in  the  30°-75°  edge  model  using  the  preceding  parameters.  Damage  values 
corresponding  to  the  contour  levels  are  also  tabulated  in  these  figures.  At  an  applied  load 
of  26.4  MPa  [Fig.  4.28(d)],  the  damage  within  the  ligament  area  w-as  still  -0.25  which  is 
less  than  the  saturation  value  of  0.53.  Since  the  observed  coalescence  loads  for  the  30°- 
75°  flaws  were  15.7/17.8  MPa  (see  Section  3.3.3),  it  was  apparent  that  the  initial  choice  of 
parameters  was  quite  inappropriate.  Specifically,  the  saturation  stress,  o^t  of  35  MPa  was 
probably  too  high  and  should  be  lowered  by  decreasing  both  A  and  £0  [see  Eqn.  (4.26)]. 
This  was  achieved  in  the  succeeding  choice  of  parameters. 

In  the  second  run,  rather  than  setting  e0  to  0.003,  the  strain  at  damage  saturation 
E<,at  was  set  to  0.003,  an  assumption  which  effectively  lowered  the  value  of  e0.  A  second 
condition  was  required  to  uniquely  define  the  two  unknown  parameters.  Since  in  the  initial 
run  it  was  observed  that  was  too  high,  in  the  second  run  this  was  prescribed  as  a  factor 
of  Ee0  as  follows: 


oui  =  1.1Ee0 


(4.29) 


This  implies  that  damage  saturation  is  reached  at  a  stress  which  is  only  107c  above  the 
stress  at  which  the  damage  initiates,  i.e.,  Ee0.  By  combining  Eqn.  (4.26)  and  (4.29),  and 
since  esa[  =  e^A,  e0  and  A  were  found  from  the  following  equations: 


o 


m: 


Eg0  1 

A  exp(l-A) 


l.lEe0 
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e  ,  =  — =  0.003 
5  A 

The  parameters  which  satisfy  the  above  equations  are:  Eo  =  0.00175  and  A  =  0.583,  with 
Dsat  =  0.34.  The  resulting  stress-strain  curve  is  shown  in  Fig.  4.29.  Note  that  is  -13 
MPa.  less  than  half  of  the  value  for  the  initial  run. 


The  damage  evolution  with  increasing  load  for  e0  =  0.00175  and  A  =  0.583  is 
shown  in  Figs.  4.30(a)-(d).  At  an  applied  load  of  1 1  MPa  (Fig.  4.33(d)),  the  region  of 
damage  with  D  =  D,,al  has  connected  the  flaw  tips.  Therefore,  according  the 
aforementioned  coalescence  criterion,  the  coalescence  load  predicted  by  £<,  =  0.00175  and 
A  =  0.583  is  ~  1 1  MPa.  The  fact  that  the  coalescence  load  for  this  set  of  parameters  is 
lower  than  the  value  computed  from  the  initial  run  is  not  surprising  given  the  smaller  value 
of  osat  in  the  second  run.  However.  1 1  MPa  is  lower  than  the  observed  coalescence  loads 
of  15.7  and  17.8  MPa.  The  assumed  o^t  was  probably  too  low;  consequently,  a  higher 
Osat  was  prescribed  for  the  succeeding  run. 


For  the  third  simulation,  Esat  was  still  set  to  0.003  as  in  the  previous  run. 
However,  osa(  was  set  to  a  value  of  15  MPa,  a  value  which  is  slightly  higher  than  13  MPa 
for  the  previous  run.  Thus,  the  parameters  for  the  third  simulation  were  found  from  the 
following  equations: 


o 


m: 


E£o 

A 

E... 


- =  15  MPa 

exp(l  -  A) 


=  ^  =  0.003 
A 


The  parameters  which  satisfy  the  above  equations  are:  Eo  =  0.0024  and  A  =  0.785,  and  the 
resulting  stress-strain  curve  is  shown  in  Fig.  4.31.  Due  to  a  slightly  higher  prescribed 
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o^.  this  new  choice  of  parameters  w  as  expected  to  predict  higher  coalescence  loads  which 
could  be  closer  to  the  observed  values. 

Figs.  4.32(a)-(d)  show  the  damage  evolution  for  to  =  0.0024  and  A  =  0.785,  with 
Dsat  =  0.19.  In  Fig.  4.32(d)  where  the  applied  load  is  13.6  MPa,  there  is  an  entire  region 
between  the  flaw  tips  where  the  damage  is  approximately  equal  to  D^.  Thus,  the  predicted 
coalescence  load  is  -13.6  MPa  and  was  thought  to  be  sufficiently  close  to  the  observed 
values  of  15.7/17.8  MPa.  Subsequent  analyses  on  other  flaw  geometries  were  performed 
using  e0  =  0.0024  and  A  =  0.785. 


Model  Validation 

Three  other  flaw  geometries,  namely  30°-45\  45°-60°  and  60°-45°  flaws,  for 
which  coalescence  was  observed  in  the  experiments  were  analyzed  using  the  strain-based 
damage  model  w  ith  t0  =  0.0024  and  A  =0.785  with  Dsal  =  0.19.  The  simulated  damage 
evolution  for  each  geometry  are  shown  in  Fig.  4.33(a)-(d),  4.34(a)-(c)  and  4.35(a)-(d). 
In  all  three  sets  of  figures,  the  last  figure  shows  a  region  between  the  flaw  tips  where  D  = 
Dsa[  such  that  the  loads  associated  with  these  contour  plots  were  chosen  as  the  coalescence 
loads  (e.g.  Fig.  4.33(d)  for  30°-45°,  Fig.  4.34(c)  for  45°-60°  and  Fig.  4.35(d)  for  60°- 
45°).  The  predicted  and  observed  coalescence  loads  for  all  four  geometries  which  were 
analyzed  are  listed  in  Table  4.2.  These  are  also  plotted  against  ligament  orientation  angle  in 
Fig.  4.36,  where  the  hollow  symbols  represent  observed  coalescence  loads  and  the  filled 
symbols  represent  predicted  coalescence  loads.  There  is  a  satisfactory  agreement  between 
the  experimental  and  computed  values  with  discrepancies  that  are  roughly  the  same  for  all 
geometries.  An  additional  calibration  cycle  would  have  resulted  in  an  even  better  agreement 
between  experimental  and  computed  coalescence  loads.  However,  the  current  results 
already  show  that  the  trends  in  coalescence  loads  as  a  function  of  flaw  geometry  were  very 
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captured  by  the  smeared-crack/strain-based  damage  criterion  proposed  in  this  study.  Any 
additional  parameter  refinement  would  just  be  an  exercise  in  curve  fitting. 

Table  4.2  Comparison  between  Experimental  and  Computed  Coalescence  Loads 


Flaw  Geometry 

Coalescence  Loads  (MPa) 

Experimental 

Computed 

30°-75° 

15.7/17.8 

13.6 

30° -45° 

24.2  /24.1 

21.4 

45°-60° 

22.4  /  23.3 

18.0 

60°-45° 

29.6  /  30.3 

24.5 

One  final  simulation  was  performed  on  30° -30°  flaws  for  which  no  coalescence  was 
observed.  Figs.  4.37(a)-(c)  show  the  damage  evolution  for  the  edge  model  of  the  30°-30° 
flaws.  At  an  applied  load  of  23.4  MPa,  damage  zones  extending  away  from  the  ligament 
had  substantially  grown  but.  in  contrast  to  the  previous  simulations,  damaged  regions  did 
not  "connect"  within  the  ligament  zone.  This  is  in  agreement  with  the  fracturing  that  was 
observed  for  the  30°-30°  flaws.  As  mentioned  in  Section  3.3.1,  one  of  the  blocks  "failed" 
at  an  applied  load  of  30  MPa  through  secondary'  crack  growth  but  these  secondary  cracks 
did  not  form  a  coalescence  crack  (see  Fig.  3.32).  Thus,  the  strain-based  damage  model 
was  also  successfully  used  to  predict  the  non-occurence  of  coalescence  between  30°-30° 
flaws. 
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Figure  4.27  Stress-strain  curv  e  for  initial  choice  of  parameters  to  simulate  coalescence, 
fo  =  0.003  and  A  =  0.25. 
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Stress-strain  curve  for  second  choice  of  parameters  to  simulate  coalescence. 
e0  =  0.00175  and  A  =  0.583. 
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Figure  4.36  Comparison  of  experimental  and  computed  coalescence  loads  for  non¬ 
overlapping  flaws. 
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Chapter  5 
Conclusions  and 

Recommendations  for  Future  Study 


An  extensive  series  of  crack  growth  experiments  using  gypsum  as  a  model  material 
for  rock  showed  fracturing  that  was  quite  different  from  that  observed  by  others  in  polymer 
materials.  In  particular,  only  wing  cracks  grew  from  pre-existing  flaws  in  polymer 
specimens  w  hile  in  gypsum,  w  ing  cracks  grew  first  and  were  follow-ed  by  secondary  crack 
growth.  More  importantly,  secondary  crack  grow'th  led  to  coalescence  between  non¬ 
overlapping  flaws  in  gypsum  specimens  w’hereas  such  coalescence  has  never  been 
observed  in  previous  experiments  on  polymers.  On  the  other  hand,  coalescence  between 
overlapping  flaws  in  gypsum  occured  through  wing  crack  growth  and  at  loads  that  were 
higher  than  the  coalescence  loads  measured  for  non-overlapping  flaws.  Thus,  analytical 
modeling  in  this  study  was  focused  on  representing  secondary  crack  growth  and  the 
associated  coalescence  between  pre-existing,  non-overlapping  flaws. 

In  the  analytical  phase,  linear  elastic  stress  and  strain  fields  after  w'ing  crack  grow  th 
were  determined  for  various  interacting  flaw  geometries.  Principal  stress  and  strain  plots 
from  these  analyses  revealed  that  secondary  crack  growth  can  be  best  modelled  by  a 
smeared  crack/damage  approach  coupled  with  a  tensile  strain-based  failure  criterion.  A 
non-softening,  tensile  strain-based  damage  model  was  then  developed  and  implemented  as 
a  material  subroutine  in  a  finite  element  program  (ABAQUS).  Material  parameters  were 
chosen  such  that  the  simulated  coalescence  load  for  one  of  the  flaw  geometries  was  close  to 
the  experimental  value.  To  validate  the  model,  the  same  set  of  parameters  was  used  to 
simulate  coalescence  between  flaws  with  other  arrangements.  Very  satisfactory  agreement 
was  achieved  between  experimental  and  simulated  coalescence  loads. 
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Although  the  damage  model  presented  in  this  study  was  derived  from  tests  on 
several  different  flaw  arrangements,  it  would  be  interesting  to  validate  this  model  on  other 
flaw  geometries  which  may  be  obtained  by  varying  the  flaw  inclination  angle,  ligament 
(rock  bridge)  inclination  angle,  flaw  length  and  ligament  length.  Other  variations  of 
specimen  geometry  are  specimens  with:  (1)  more  than  tw'o  flaws,  (2)  sub-parallel  or  even 
non-parallel  flaws,  and  (3)  closed  flaws.  The  results  of  these  additional  fracture  tests  could 
then  be  compared  to  results  of  numerical  simulations.  It  would  also  be  interesting  to  see 
whether  the  fracturing  of  double-edge  notched  concrete  beams  described  in  Chapter  1  of 
this  report  could  be  simulated  using  the  tensile  strain-based  damage  model.  In  fact,  since 
the  model  is  based  on  stresses  and  strains,  it  could  be  applied  to  any  geometry  as  long  as 
the  material  remains  elastic  prior  to  fracture,  and  the  fracturing  of  the  material  is  adequately 
characterized  by  a  strain-based  failure  criterion. 

Another  obvious  extension  of  this  work  would  be  to  perform  the  crack 
growth/coalescence  experiments  on  other  types  of  rock  and,  eventually,  different  brittle 
materials.  Since  secondary  crack  growth  was  seen  in  other  rocks  as  reported  in  the 
literature,  it  is  highly  possible  that  coalescence  would  also  occur  in  these  other  rock  types 
as  it  did  in  gypsum.  Different  rock  types  have  different  microstructures,  and  it  would  be 
interesting  to  see  how  the  coalescence  behavior  is  affected  by  these  differences. 

In  this  study,  the  stress-strain  distribution  around  the  flaws  in  the  test  specimens 
was  only  established  through  numerical  simulation  of  the  experiments.  The  non¬ 
homogeneity  and  strong  gradients  of  the  strain  field  around  the  flaws  made  point-strain 
measurement  techniques  such  as  using  strain  gages  infeasible.  However,  techniques  are 
currently  being  developed  for  measuring  "whole  strain  fields"  in  non-photoelastic  materials 
such  as  rock.  Examples  of  these  are  holographic  and  speckle  interferometry.  The 
application  of  these  methods  to  the  crack  growth/coalescence  experiments  described  in  this 
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report  would  not  only  verify  the  numerically  derived  stress  and  strain  fields  but  would  also 
confirm  or  refute  the  hypothesis  that  fracturing  occurs  in  regions  w-ith  highest  tensile 
strains. 
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